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ABSTRACT 


The  purpose  of  this  thesis  was  to  evaluate  the  accuracy  of  the  existing  Mini- 
Ranger  network  in  Monterey  Bay  and  to  suggest  means  by  which  the  network 
can  be  improved.  This  network  consists  of  six  stations  located  around  the  Bay 
and  installed  by  the  Monterey  Bay  Aquarium  Institute  (MBARl)  and  the  Naval 
Postgraduate  School  (NPS). 

In  order  to  undertake  this  task,  data  were  made  available  by  MBARl  which 
were  collected  on  a  cruise  made  by  their  vessel  "Pt.  Lobos".  Additional  data  were 
gathered  on  a  second  cruise  by  the  vessel  "Pt.Sur"  of  the  NPS. 

The  data  analysis  indicated  network  problems.  An  effort  was  made  to  iden¬ 
tify  these  problems  and  to  compute  various  correctors.  In  addition,  an  equation 
has  been  derhed  which  enr.b'es  use  of  Mini-Ranger  data  collected  when  signal 
strengths  are  low. 

The  estimation  of  the  accuracies  obtained  from  the  network  through  the  \  ar- 
ious  tests  applied,  and  the  conclusions  drawn,  can  be  used  as  a  guide  to  future 
users  of  the  system. 
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I.  INTRODUCTION 


In  1988,  the  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  and  the 
Naval  Postgraduate  School  (NPS)  began  work  to  create  a  permanent  network  of 
Mini-Ranger  stations  around  Monterey  Bay.  The  net  was  to  provide  precise  po¬ 
sitioning  for  research  vessels  and  other  scientific  operations.  Initially  MBARI 
was  seeking  coverage  of  as  much  of  the  bay  as  possible,  especially  the  Monterey 
canyon  feature,  with  absolute  horizontal  position  acuracies  on  the  order  of  a  few 
meters.  It  was  understood  this  would  be  difficult  to  achieve,  gi\en  the  equip¬ 
ment's  specified  range  (40  nm).  This  study  was  undertaken  to  determine  what 
levels  of  accuracy  and  spacial  coverage  were  being  achieved,  as  the  net  was  con¬ 
figured  in  early  1989.  and  also  to  suggest  possible  impro\ements. 

Positioning  data  from  two  cruises  were  examined,  one  from  a  "Point  Lobos" 
cruise  on  31  MAR  89,  and  a  second  from  "Point  Sur"  cruise  on  22  SEP  89  (the 
"Point  Lobos"  and  "Point  Sur"  are  vessels  used  by  MBARI  and  NPS.  respec- 
ti\ely).  Significant  errors  were  found  in  a  high  percentage  of  the  position  fixes 
from  the  "point  Lobos"  data  set.  It  was  apparent  that  at  least  some  of  the  prob¬ 
lem  arose  from  the  use  of  weak  signals  (low  signal  strength  readings)  and  incom¬ 
pletely  calibrated  equipment.  Consequently,  additional  measurements  were  made 
onshore  o\er  known  length  baselines  in  order  to  quantify  the  systematic  and 
random  error  components  of  the  range  measurement  under  \'arious  signal 
strength  conditions. iThese  findings  stimulate  the  direction  of  the  rest  of  the 
study: 

•  What  were  the  main  sources  of  errors? 

•  Is  there  any  postprocessing  that  can  correct  the  data  and  derive  better  posi¬ 
tions? 

•  What  is  a  realistic  estimate  of  the  accuracy  achieved? 

•  Are  there  any  chances  for  impro\ement? 

1  It  should  he  noted  that  the  "Point  Lobos"  Master  station  failed  before  the  baseline  testing 
was  run.  I  hus  it  was  not  possible  to  confirm  the  error  estimates  derived  from  analysis  of  the  cruise 
data  with  independent  baseline  measurements. 
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II.  DESCRIPTION 


A.  MINI-RANGER  SYSTEM 

The  Motorola  Mini-Ranger  system  being  used  by  MBARI  and  NFS  is  the 
Falcon  484.  It  is  widely  used  for  precise  positioning  in  a  variety  of  nearshore 
marine  applications  like  hydrographic  surveying,  dredging,  oceanographic  data 
collection,  etc.  Numerous  excellent  descriptions  of  the  system  are  avaliable  in  the 
literature. 

The  system  operates  on  the  principle  of  a  pulse  coded  transponder.  The 
Master  station  (also  refered  as  to  as  the  Receiver-Transmitter)  is  located  on  the 
ship  and  ranges  to  a  maximum  of  four  Reference  stations  ashore.  The  maximum 
update  rate  is  once  per  second.  The  Master  and  Reference  stations  transmit  on 
different  frequencies  in  the  microwave  band.  5570  and  5480  Mhz,  respecii\’ely. 
Both  range  and  signal  strength  data  are  output  for  each  Reference  station  inter¬ 
rogated. 

Signal  strengths  are  computed  at  the  Master  as  a  re]ati\e  indication  of  the 
power  level  received  from  the  Reference  station.  Values  can  vary  from  about  5 
to  99.  The  manufacturer  warns  that  any  ranges  with  signal  strength  below  13  are 
"weak"  signals,  subject  to  random  errors  exceeding  the  2  meters  (1  cr  )  specifica¬ 
tion.  Signal  strength  obx  iously  falls  off  with  increasing  range,  but  can  also  drop 
due  to  any  degradation  of  propagation  conditions  between  the  ship  and  shore 
station.  Hence,  weak  signals  can  and  do  occur  at  any  range.  Monitoring  of  the 
received  signal  strength  is  critical  to  assessing  the  accuracy  of  all  range  data. 

The  range  is  computed  at  the  Master  station  based  on  the  measured  round- 
trip  travel  time  for  the  pulse  and  a  preset  refractive  index  for  microwave  propa¬ 
gation  in  the  atmosphere.  Unless  changed  by  the  user,  the  recei\cr  assumes  a 
value  of  N  =(n-l)xl0^  =  320  as  an  indicator  for  the  refractive  index  (n).  which  is 
the  value  used  throughout  this  study.  Output  ranges  also  can  be  corrected  for  any 
constant  (systematic)  error  previously  entered  by  the  user.  Because  this  constant 
error  includes  time  delays  at  both  ends  of  the  line,  it  is  essential  to  haw  a  specific 


calibration  of  each  Master/ Reference  station  pair.  The  manufacturer  notes  that 
without  calibration,  systematic  range  errors  of  up  to  ten  meters  are  possible. 

It  is  possible  to  ha'.e  the  receiver  correct  ranges  from  "slope  to  horizontal" 
distances  before  output,  based  upon  user  provided  elevations  of  the  Master  and 
Reference  stations.  During  this  study,  both  the  MBARI  and  NFS  equipment 
were  configured  to  output  "slope"  ranges  only,  choosing  to  handle  this  correction 
as  part  of  the  position  computation  algorithm. 

The  Reference  stations  are  supplied  with  electric  current  of  22  to  32  Volts 
DC  from  a  pair  of  common  car  batteries,  or  permanently  from  an  AC  supply 
through  a  con\erter.  A  possible  insufficient  current  supply  can  cause  errors  in  the 
displayed  distance  in  the  receiver  [Mini-Ranger  Operation  Manual,  p.11-1 1]. 

The  standard  system  has  a  maximum  range  of  37  km  (20  n.m),  with  Refer¬ 
ence  station  antennas  of  13  db.  The  maximum  range  can  be  extented  to  75  km 
(40  n.m)  by  raising  the  antenna  gain  of  the  Reference  stations  to  19  db,  if  the 
station's  elet  ation  is  sufficient  (due  to  the  dip  of  the  horizon).  The  master  station, 
which  is  on  the  t'essel.  has  an  omni  directional  antenna  of  6  db.  The  discussion 
and  results  in  this  thesis  are  all  given  using  extended  (19  db)  gain  antennas  on  the 
shore  stations. 

Another  significant  factor  for  a  Reference  station  is  the  dirccth  iiy  of  its  an¬ 
tenna.  The  gi\  en  maximum  range  of  the  system  corresponds  to  the  main  lobe  di¬ 
rection  and  drops  as  one  mo\  es  to  the  sides.  Figure  1  shows  the  pattern  of  a  19 
db  gain  antenna  of  a  shore  station.  We  can  see  that  the  range  drops  under  30 
km  at  40^  off  the  center  direction.  This  is  significant  for  a  station  that  uses  its 
maximum  range  ability  and  thus  makes  the  choise  of  antenna  directit  ity  eritical. 

B.  MONTEREY  BAY  NETWORK  DESCRIPTION 

The  Motorola  Mini-Ranger  network  in  Monterey  Bay  consists  of  6  stations 
with  high  gain  antennas  (19  db),  installed  permanently  on  known  positions 
around  the  Bay.  The  installation  was  done  by  both  MBARI  and  NFS  serxing 
primarily  the  Biological  and  Oceanographic  research  of  their  xessels  "Foint 
Lobos"  and  "Foint  Sur"  respectixely.  The  positions  at  xxhich  the  stations  are  in- 
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Figure  1.  19-db  ANTENNA  PATTERNS. 
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Table  I.  PRELIMINARY  POSITIONS  OF  THE  MINI-RANGER 
STATIONS 


Name 

Latitude(N) 
Longitude(W) 
deg  min  sec 

UTM  coordi¬ 
nates 

Ele¬ 

vation 

(m) 

T  revor 

36  58  28.785 

122  02  31.536 

585260.161 

4092490.284 

.999690 

51.7 

Wats 

36  54  26.895 

121  50  39.819 

602947.672 

4085231.399 

.999731 

23.1 

Packard 

36  49  44.407 

121  46  04.912 

609863.128 

4076611.345 

.999749 

33.9 

Hays 

36  38  33.826 

121  47  45.895 

607621.289 

4055915.264 

.999743 

137.2 

Doppler 

36  36  05.601 

121  52  37.479 

600434.009 

4051260.014 

.999724 

10.0 

Hank 

36  36  25.381 

121  55  08.722 

596669.490 

4051826.411 

.999715 

134.5 

stalled  are  not  ideal  due  both  to  the  permanent  power  supply  problem  and  to  the 
fact  that  some  of  the  ideal  station  positions  are  on  pri\  ate  properties. 

On  board  the  "Point  Lobos"  the  Mini-Ranger  omni-directional  antenna  was 
installed  on  the  mast  9.1  meters  abo\e  the  water  le\el.  Thus  the  maximum 
achie\  ablc  range  due  to  the  line  of  sight  effect  was  increased. 

The  shores  of  the  Bay  and  the  nearshore  areas  arc  generally  low  sand  dunes 
so  some  stations  don't  ha\e  the  proper  altitude  to  reach  their  maximum  range. 

Table  1  gi\es  the  preliminary  positions  [Schnebcle,  1989]  of  the  stations,  the 
UTM  grid  coordinates  and  their  scale  factor.  The  horizontal  datum  is  the 
American  Datum  of  1983  (NAD  83).  The  UTM  projection  used  is  for  Zone  10 
(central  meridian  123'  W.)  and  the  scale  factor  0.9996.  Any  errors  in  the  posi- 


5 


lions  of  these  stations  will  propagate  into  the  computed  vessel  position  by  an 
amount  that  depends  on  the  geometry.  If  the  errors  are  smaller  than  several 
decimeters  then  there  is  no  significant  effect.  Revised  positions  are  shown  in  Ta¬ 
ble  3,  section  IV-C. 

The  stations  are  not  ideally  positioned  and  so  may  be  moved  to  different  po¬ 
sitions  in  the  future  to  allow  belter  coverage  of  the  areas  of  interest. 

Station  Doppler  is  not  yet  permanently  installed  for  safety  reasons. 
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III.  THEORETICAL  APPROACH  TO  THE  PROBLEM 
A.  THE  RANGING  PROBLEM 

The  technique  used  for  computing  the  ship's  position  using  ranges  from 
known  reference  stations  will  be  derived  and  explained  in  the  beginning  so  the 
reader  will  be  able  to  follow  and  understand  the  data  manipulation  and  results 
of  this  study.  It's  well  known  that  one  way  to  handle  the  computations  invoh  ing 
ranging  data  is  to  work  on  some  appropriate  projection  (e.g  the  UTM  projection). 
There  are  several  techniques  that  can  be  used  for  the  position  derivation: 

•  At  eraging  the  coordinates  resulting  from  all  possible  combinations  of  meas¬ 
ured  ranges. 

•  Position  calculation  using  the  two  most  accurate  ranges  (using  of  the  rest  as 
a  check). 

•  Graphical  methods. 

In  this  thesis  we  use  the  variation  of  coordinates  method  which  is  described 
in  Cross  [1981].  For  the  position  derivation,  at  least  2  ranges  are  needed  (two 
unknown  parameters  in  the  plane  solution  X.Y).  With  more  ranges  one  can  esti¬ 
mate  the  accuracy  of  the  resulting  position  and  also  make  statistical  decisions 
about  the  quality  of  the  data. 

Figure  2  illustrates  the  concept  of  ranging  from  the  i-th  station.  Consider  a 
station  S,  in  a  coordinate  system  XY.  The  observed  range,  reduced  to  the  plane, 
is  OR,  (not  shown).  Consider  an  assumed  position  (AP).  The  distance  between 
the  reference  station  and  the  assumed  position  is  CR,.  The  calculation  of  the 
values  dX.dY.  which  represent  the  corrections  that  ha\e  to  be  applied  on  the  AP, 
will  giN  e  the  desired  true  position  (TP)  after  some  iterations. 

As  is  ob\  ious  from  the  figure: 

^  AP~  ^’si 

The  key  equation  is  the  one  for  the  calculated  range: 
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Figure  2.  VARIATION  OF  COORDINATES  METHOD 


In  order  to  find  the  transformation  from  range  coordinates  to  position  coor¬ 
dinates  it  is  necessary  to  find  the  derivative: 


ICR^dCRi  =  2{X^^pdX  -f  2( 


dCRi  =  ^^dX  -t-  NdY 

where 
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CRi 


dY 


The  estimate  of  the  true  range  then  is  equal  to  the  observed  range  plus  an 
error  which  can  be  positive  or  negative,  called  the  residual  V,.  It  is  also  equal  to 
the  calculated  range  plus  the  increment  dCR,. 

The  true  range  then  is  OR^  +  V,  where  F,  is  the  residual.  But  it's  also 
CR,  4-  dCR,.  Solving  for  dCR^,  we  take: 

dCRi  =  (OR  -  CR)i  +  I'i 


and  finally; 


MidX  +  N,dY  =  (OR  -  C/?),  +  1 

For  the  n  ranges  problem,  in  matrix  form,  the  equation  may  be  slated  as  fol¬ 
lows; 


AX  =B+V 


Where  matrices  A  ,B  and  X  respcctiscly  are: 


A  = 


.W, 

M2 


^2 


M. 


(OR  -  cy?), 
(OR  -  CR)2 

(OR  -  CR)„ 


9 


The  standard  solution  for  this  formulation  using  the  least  squares  method 
(Bomford,  1980)  gives: 

a;,  =  {A  ^IVAf^A  ^WB 
l\=^AX^-B 

Where  the  subscript  a  indicates  that  the  value  is  approximate.  The  best  sol¬ 
ution  is  obtained  by  iterating  until  X  converges  to  less  than  some  limit. 

W  is  the  weight  matrix.  The  standard  daviation  tr,  of  a  range  measurement 
is  an  e.xpresion  of  its  accuracy.  It  is  standard  to  weight  the  range  measurements 
using  -"Y  which  ensures  achieving  the  minimum  variance  solution  [Hamilton. 
1964]. 

Assuming  that  there  is  no  correlation  between  measurements,  the  weight 
matrix  becomes  diagonal  as  follows: 

»'=  “ 

0  0 

Motorola  estimates  the  standard  deviation  (1  cr  )  of  a  single  range  measure¬ 
ment  to  be  2  meters.  The  individual  user  has  to  set  his  own  estimated  (  o  )  for 
his  equipment  (see  section  IV-H). 

The  unbiased  estimate  of  the  variance-covariance  matrix  of  the  adjusted  po¬ 
sition  is  given  by 

This  matrix  includes  the  variances  and  the  covariances  in  the  coordinate  sys¬ 
tem  we  use  (  X,Y  coordinates)  and  appears  as  follows: 
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(18) 


Assuming  that  a,  =  2m,  then  the  combination  of  two  ranges  can  result  in  po¬ 
sitional  error  of  7.8  m  (drms)  in  the  limiting  case  of  30"  intersection  angle,  Drms 
is  derived  from  the  unbiased  estimate  of  the  variance  covariance  matrix  as  fol¬ 
lows: 

\drms  =  +  cj) 

This  error  of  7.8  m  can  be  larger  when  the  system  is  not  calibrated.  For  un¬ 
calibrated  systems  the  error  of  a  single  measurement  can  exceed  10  m  [Motorola 
Mini-Ranger  manual]. 

The  quantity  cr*  is  the  a  posteriori  variance  of  unit  weight.  It  can  be  computed 
as  follows: 


2  (I'Jn'i;) 

do/ 

Where  dof  is  the  degree  of  freedom,  or  n-2  when  solving  for  X.^'  position  co¬ 
ordinates  using  n  ranges.  Ideally.  al  should  equal  unity. 

According  to  Cross  [1981],  when  the  degrees  of  freedom  are  small  it  is  dan¬ 
gerous  to  use  a  single  determination  of  to  judge  the  quality  of  the  position.  It 
is  better  to  use  an  average  value  of  <7^  from  the  whole  data  set. 

But  the  big  question  is  "how  we  can  say  that  the  value  of  is  out  of  bounds? 
For  this  question  a  chi  square  one  or  two  tailed  test  is  the  best  indication,  even 
though  there  are  situations  where  it  is  not  sufficient  [Uotila.  1975].  If  the  test 
fails,  this  may  be  an  indication  that  one  or  more  of  the  following  error  sources 
exist  that  cause  the  large  s  alue  of  g^. 

•  Error  in  the  mathematical  model  involved  in  the  solution  (Scale  factors. 

Refractive  Index  used  etc.) 

•  Computational  errors. 
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•  High  correlations  between  observations  or  very  poor  geometry  of  the  situ¬ 
ation. 

•  Influence  of  the  omitted  higher  order  terms  in  Taylor  series  expansion  which 
can  be  critical  for  higher  order  accuracy  systems  (primarily  when  the  Master 
station  is  far  away  from  to  the  Reference  station). 

•  Incorrect  variance-covariance  matrix  of  the  observed  quantities.  It  means 
that  for  a  non  well  calibrated  system,  large  values  of  are  expected. 

•  Blunders  in  observations  arising  from  various  factors  in  the  measurement 
process  (Multipath,  unstable  readings,  etc). 

So  if  the  statistical  test  used  for  shows  that  something  is  wrong  an  in\  esti- 
gation  is  needed  for  the  detection  of  the  possible  error  source  and  its  elimination 
from  the  data. 

B.  PROPAGATION  EFFECTS 

The  line  of  sight  limitation  is  one  of  the  critical  factors  in  the  network  instal¬ 
lation.  The  radius  of  cur\  ature  of  the  signal  path  is  about  4  times  greater  than  the 
Earth's  curvature,  so  the  heights  of  the  Reference  stations  and  the  height  of  the 
Master  station  antenna  govern  the  maximum  range  of  the  system,  according  to 
the  formula,  [Motorola  Mini-Ranger  Manual,  1981]: 

^  =  4.04(^7;7  +  nM 

Where  h,  and  h„  are  the  heights  in  the  Reference  and  Master  stations  respectively 
and  d  the  maximum  range  in  kilometers. 

The  maximum  range  of  the  system  can  be  reduced  due  to  the  atmospheric 
conditions.  Attenuation  of  the  signal  due  to  rain  is  given  by  [Casey,  1982]: 

A  =  ar^ 

Where  A  is  the  attenuation  of  the  signal  in  db,  r  is  the  rainfall  rate  in  mm  hr  and 
a.b  are  frequency  functions.  For  the  Miniranger  mean  frequency  (5.5  GHz),  their 
numerical \alues  are  a  =  1.48  x  10'^  and  b=  1.1469. 

Another  formula  is  gi\  en  for  fog.  based  on  the  water  concentration  per  cubic 
meter: 
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A  =  4.87xlO~^M/db/km. 


Where  M  is  the  water  content  in  grjm^  and  f  is  the  frequency  in  Ghz. 

Table  2  gives  the  attenuation  of  the  signal  for  rain  and  fog,  at  5.5  Ghz  which 
is  an  average  operating  frequency  for  the  Mini-ranger.  In  terms  of  distance,  we 
have  8.0  km  range  reduction  at  40  km  for  each  db  drop.  So,  for  example,  in  a 
distance  of  40  km  with  heavy  rain  conditions  we  have  1.2  db  drop  in  the  signal 
(1 1.6  km  reduction  in  maximum  range).  This  drop  in  signal  power  can  introduce 
unexpected  range  errors  when  signal  strength  drops  below  the  Critical  Strength 
Threshold  (CST). 


Table  2.  ATTENUATION  OF  THE  SIGNAL  AT  5.5  MHZ 


Precipitation 
Rate  (mm  hr) 

Attenuation 
(db  km). 

Fog  (gr  /»■') 

Attenuation 
(db  km) 

2 

0.0015 

1.0 

0.0147 

4  (drizzle) 

0.0072 

1.2 

0.0177 

6 

0.01 16 

1.4 

0.0206 

8 

0.0161 

1.6 

0.0236 

10 

0.0208 

1.8 

0.0265 

12 

0.0256 

2.0 

0.0295 

14 

0.0305 

2.2 

0.0324 

16 (heavy) 

0.0356 

2.4 

0.0354 

18 

0.0407 

2.6 

0.0383 

20 

0.0460 

2.8 

0.0413 

A  possible  problem  that  can  appear  in  the  survey  area  is  the  multipath  phe¬ 
nomenon  of  the  signal.  It  can  cause  constructive  or  destructive  inteferencc  at  the 
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receiver.  There  exist  two  types  of  multipath.  The  first,  vertical  multipath  is 
caused  by  reflection  of  the  signal  on  the  sea  surface.  The  destructive  interference 
at  the  receiver  due  to  vertical  multipath  is  called  "range  hole".  The  second  one 
is  the  horizontal  multipath  due  to  reflection  of  the  signal  on  walls  or  other  sur¬ 
faces.  The  destructive  interference  of  the  horizontal  multipath  is  seen  as  "unstable 
readings". 

Range  holes  are  very  important  in  networks  that  operate  at  line  of  sight  and 
can  be  a  major  source  of  positional  problems.  The  best  way  to  avoid  such  a 
problem  is  to  put  an  additional  Mini-Ranger  Master  station  at  a  different  height 
on  the  vessel. 

The  accuracy  of  the  range  measurements  is  affected  by  the  weather  conditions 
(temperature,  humidity,  pressure).  A  fixed  value  of  the  refractis  e  index  (N)  is  in¬ 
put  into  the  system  prior  to  a  survey  and  kept  constant  throughout  the  sur\ey. 
The  default  refractixe  index  (N=320)  corresponds  to  atmospheric  conditions  of 
20^C,  50®  0  and  1013  mb  respectively. 

Extreme  atmospheric  conditions,  like  high  temperature  along  with  either  high 
relatixe  humidity  or  dry  atmosphere,  can  give  rise  to  large  scale  errors  in  the 
measured  distances.  For  this  reason  further  investigation  was  done  to  estimate 
the  magnitute  of  the  atmospherically  induced  errors  (see  section  IV-E). 

C.  LOW  SIGNAL  STRENGTH 

Sexeral  investigators  [e.g  Casey,  1982]  have  shown  that  range  errors,  both 
bias  and  standard  deviation,  tend  to  increase  once  the  received  signal  strength 
drops  below  some  critical  level.  This  level,  usually  refered  to  as  the  Critical  Signal 
Threshold  (CST).  is  defined  by  the  users  accuracy  requirements.  Typical  choices 
are  between  signal  strengths  of  8  to  18.  The  Motorola  Falcon  equipment  itself 
flags  signal  strengths  less  than  or  equal  to  13  as  potentially  suspect. 

The  signal  strength  from  the  shore  station  is  also  a  factor  that  can  introduce 
large  errors  in  the  range  measurements,  but  it  is  always  output  bv-  the  svstem  and 
can  be  used  in  subsequent  data  processing.  When  signal  strength  is  above  the 
CST  limit,  the  errors  have  a  Gaussian  distribution,  (mean  v  alue  zero  and  rela¬ 
tively  small  variance).  But  when  the  signal  strength  is  below  the  CST  limit  tlie 
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errors  follow  a  logarithmic  curve.  That  means  the  range  error  increases  logarith¬ 
mically  as  one  recci\es  smaller  values  of  signal  strength  [Casey,  1982]. 

Figure  3  derixed  by  Casey,  gives  the  CST  at  SS  =  8  .  For  particular  applica¬ 
tions  each  user  can  set  his  own  estimated  value  for  CST  in  order  to  cither  reject 
the  low  signal  strength  meausurements.  or  to  increase  the  standaid  dexiaiion  of 
the  measurement. 

The  SS  of  each  signal  depends  strongly  on  the  distance  from  the  rcrcrcncc 
station,  (it  drops  with  the  inxerse  square  of  the  distance).  It  can  also  \  ary  due  not 
only  to  multipath  phenomenon,  (as  we  explained  earlier),  but  also  due  to  the  an¬ 
gle  to  the  main  lobe  direction  and  the  atmospheric  conditions. 
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IV.  TESTS  ON  THE  MINI-RANGER  NETWORK  IN  MONTEREY  BAY 

A.  MINI-RANGER  CRUISE  DATA 

On  31  March  1989  the  Vessel  Point  Lobos,  from  MBARI,  sailed  for  a  trip  in 
Monterey  Bay  to  make  trials  and  calibration  on  the  Mini-Ranger  stations.  The 
Master  station  was  installed  on  the  mast  at  height  10  m  above  the  water.  During 
this  trip  useful  data  were  collected  and  position  estimates  calculated.  This  data 
was  processed  in  a  preliminary  fashion  in  order  to  detect  possible  error  sources. 

Figure  4  shows  the  routes  of  the  vessel  followed.  A  division  of  the  Bay  into 
smaller  areas  was  used,  (as  shown  on  the  figure),  in  order  to  subdivide  the  data 
for  processing  and  error  anabasis  purposes.  The  division  was  based  on  the  con¬ 
cept  of  nearshore  and  open  Ln  y  regions. 

Another  trip  was  subsequently  done  for  NPS  research  purposes  on  22  Sep¬ 
tember  1989,  and  some  more  data  collected  as  shown  on  the  same  figure.  The 
master  station  for  this  cruise  was  different  and  the  data  w  collected  under  differ¬ 
ent  conditions  of  atmosphere,  time  and  instrumentation.  It  was  considered  possi¬ 
ble  that  the  positional  results  from  these  two  data  sets  may  vary  in  accuracy. 

First  we  examined  the  earlier  data  set  of  MBARI.  which  already  had  the 
positions  and  their  statistics  calculated.  We  were  informed  by  MBARI  that  for 
the  position  derivation  they  used  the  method  of  v  ariation  of  coordinates.  No  de¬ 
tails  on  the  specific  code  in  this  software  were  avaliable. 

The  first  step  was  a  general  view  on  the  residuals  and  the  "standard  devi¬ 
ation"  provided  in  MBARl's  data  of  the  whole  trip.  This  showed  that  something 
did  not  work  well  when  the  vessel  was  nearshore  ,  especially  in  the  SW  portion 
of  the  Bay.  Figure  5  shows  the  areas  of  the  greater  and  more  variable  standard 
dev  iation  of  the  data  set.  In  general  as  shown  in  Figure  5.  areas  E  and  G  have 
fewer  problems.  It  was  not  obvious  which  station  or  stations  combination  was 
responsible  for  the  big  residuals. 

Another  problem  was  the  variation  of  signal  strength  and  its  unexpectedly 
high  values.  In  many  portions  of  the  data  set,  the  signal  strength  was  much  higher 
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Figure  4.  SHIP  ROUTES  FOR  THE  TWO  DATA  SETS/AREA  DIMSION 
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than  expected  for  such  a  range.  The  appearance  of  signal  strengths  greater  than 
70  in  distances  more  than  10000  meters  is  unusual.  This  is  a  strong  indication  of 
possible  system  malfunction. 

Because  the  data  set  did  not  show  the  expected  accuracy,  the  following  series 
of  tests  and  checks  were  performed: 

•  A  visit  to  the  stations  for  inspection  of  directional  limitations,  antenna  in¬ 
stallation  problems  and  other  possible  affecting  factors. 

•  A  check  of  the  station  positions  with  geodetic  survey  techniques,  (triangu¬ 
lation,  traverse,  or  differential  GPS),  and  a  least  squares  adjustment  of  the 
survey  data  with  fixed  geodetic  control  points  around  the  Bay. 

•  Check  of  the  software  used  for  the  position  derivation  using  the  variation  of 
coordinates  algorithm  with  accurate  grid  scale  and  slope  distance  correction 
techniques. 

•  A  theoretical  view  of  the  geometrical  accuracies  with  the  use  of  various  sta¬ 
tion  combinations  in  the  Bay  area. 

•  A  check  on  the  effects  of  extreme  meteorological  conditions  to  check  for  the 
possibility  of  introducing  range  errors  in  excess  of  1  meter. 

•  A  range  hole  prediction  for  the  network. 

•  .An  investigation  into  the  possibility  of  correcting  those  ranges  with  signal 
strengths  below  the  CST. 

•  A  calibration  test  for  determining  biased  or  linearly  varying  errors  and  their 
variability  with  signal  strength  above  the  CST. 

B.  INSTALLATIO.N  CHECK  AND  ANTENNA  DIRECTIVITIES 

On  20th  of  July,  a  visit  to  the  stations  TREVOR.  W.ATS  and  P.ACK.ARD 
was  done  for  visual  examination,  check  on  the  antenna  directi\  ity  ,  possible  sig¬ 
nal  blocking  from  nearby  obstructions  and  other  possible  affecting  factors.  The 
stations  HAYS  and  HANK  were  visited  on  the  1st  of  August  and  25th  of  Sep¬ 
tember  respectively.  The  alignments  of  the  antennas  were  determined  to  approx¬ 
imately  2'-  using  a  small  portable  magnetic  compass.  All  the  magnetic  bearings 
are  gi\  en  as  true  bearings.  A  short  description  is  giv  en  below  for  each  station  and 
a  map  is  presented  (Figure  6)  for  clarification: 

•  TRE\’OR:  This  station  is  installed  on  the  roof  of  a  house  in  Santa  Cruz.  The 
antenna  is  attached  to  the  chimney  of  the  house  and  is  oriented  towards 
175''.  The  signal  is  blocked  by  a  rocky  hill  from  215  to  approximately  100' 


Figure  6.  OBSERVED  ANTENNA  DIRECTIVITIES  OF  MINI-RANGER 
NETWORK  IN  MONTEREY  BAY 
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(clockwise).  There  is  no  clear  view  of  the  Northern  near  shore  Bay  region  but 
this  area  is  of  little  interest. 

•  PACKARD;  Installed  on  the  roof  of  a  wood  craft  building,  hooked  on  a  thin 
mast,  on  private  property  between  Watsonville  and  Moss  Landing.  Its  an¬ 
tenna  points  to  250®  and  may  be  blocked  by  the  Moss  Landing  electric  plant 
stacks  between  202°  and  205°  There  were  no  other  problems  observed  for  this 
station. 

•  WATS:  Installed  on  the  roof  of  a  commercial  building  near  the  beach  at 
Watsonville.  Its  antenna  points  at  215°  .The  entire  Bay  area  is  unobstructed. 
The  antenna  is  approximately  2  feet  from  the  vertical  wall  of  the  building. 

•  HAYS;  Installed  on  the  roof  of  Silas  B.  Hays  hospital  at  Fort  Ord  with  an 
uninterupted  Bay  view  and  directivity  of  its  antenna  to  295°  .  The  perfect 
N'iew  along  with  the  height  of  the  station  makes  it  useful  in  some  portions 
SW  of  Point  Pinos  .  The  signal  is  blocked  from  the  higher  Pacific  Grove  hills 
from  253°  back  to  the  South. 

•  HANK:  Installed  on  the  roof  of  a  house  in  Pacific  Gro\e.  While  it  has  good 
height  the  signal  cannot  reach  the  nearshore  portion  from  Point  Pinos  to 
Monterey  harbour  due  to  high  trees  existing  downhill.  The  directi\  iiy  of  its 
antenna  is  345°  and  the  signal  may  been  blocked  by  high  trees  from  005°  to 
035°.  The  station  is  offset  from  the  position  specified  by  the  preliminary  co¬ 
ordinates.  This  particular  correction  is  given  in  the  discussion  in  section  C 
below. 

•  DOPPLER;  This  station  was  installed  near  NPS  Beach  Lab.  but  due  to 
safety  reasons  was  subsequenly  removed. 

It  was  not  possible  to  check  the  electric  current  supply  of  the  stations,  but  a 
voltage  stabilizer  can  help  the  system  to  be  kept  healthier  and  unstable  outputs 
can  be  a\oided  in  the  future. 

C.  POSITION  CHECK 

Geographic  coordinates  for  the  station  sites  had  been  determined  on  the 
North  American  Datum  of  1983  (NAD  83)  by  geodetic  survey  techniques,  with 
tics  to  existing  horizontal  control  points  published  by  the  National  Geodetic  Sur¬ 
vey  [Schnebele,  1989].  Conventional  ground  survey  and  Global  Positioning  Sys¬ 
tem  (GPS)  techniques  were  both  used.  The  sur\ey  data  had  been  analyzed  with 
a  three-dimcntional  least-squares  adjustment  software  package  (GEOLAB,  by 
GEOSur^■.  Inc.).  Absolute  horizontal  position  accuracies  estimated  by  the  ad¬ 
justment  were  better  than  8  cm  (circular  error,  95'’ o  confidence  le\  el).  Additional 
work  to  confirm  the  offset  positions  of  the  Mini-Ranger  units  at  P.ACKARD  and 


WATS,  has  yet  to  be  completed,  but  is  not  expected  to  shift  their  preliminary 
positions  by  more  than  a  few  centimeters  at  most. 

These  station  coordinates  were  converted  to  Universal  Transverse  Mercator 
grid  coordinates  (UTM  Zone  10)  for  use  in  the  vessel  position  computation  algo¬ 
rithms.  The  UTM  system  has  been  adopted  as  the  standard  grid  used  by  MBARl 
and  NFS  with  this  network.  The  specific  Transverse  Mercator  Projection  algo¬ 
rithm  [Floyd,  1985]  is  accurate  to  within  1  cm  for  the  entire  zone. 

As  described  earlier,  station  HANK  was  offset  from  the  preliminary  position. 
The  offset  was  measured  by  a  simple  magnetic  compass  for  azimuth  and  a  tape 
for  distance  and  should  ha\  e  errors  no  greater  than  5  cm.  The  measured  bearing 
was  066'’  true  and  the  distance  41  inches.  The  two  positions  were  at  the  same 
height,  so  the  corrections  in  UTM  are:  dx  =  0.952,  dy  =  0.424,  dz  =  0.00  in  meters. 

The  results  are  sumarized  and  the  revised  station  positions  provided  in  the 
Table  3.  With  the  exception  of  station  HANK,  the  revised  positions  of  the 
stations  differ  by  less  than  10  cm  from  the  preliminary  positions.  It  is  concluded 
that  the  station  positions  are  not  the  source  of  the  obsert  ed  inaccuracies. 

D.  ESTIMATION  AND  COMPARISON  OF  DERIV  ED  POSITIONS 
I .  Development  of  position  derivation  algorithm 

The  measured  ranges  were  used  to  recompute  the  ship's  positions  by  the 
method  of  variation  of  coordinates  [Cross,  1981].  These  recomputed  positions 
were  then  compared  with  those  provided  by  MBARl.  The  comparisons  showed 
that  there  is  a  discrepancy  between  the  given  positions  and  the  recomputed  posi¬ 
tions  on  the  first  data  set.  which  in  some  cases  were  significant.  Because  of  these 
discrepancies,  we  ignored  all  the  given  MBARl  positions  and  recomputed  all  their 
data. 

For  this  purpose,  a  program  was  written  in  Turbo  Basic.  This  program 
uses  the  variation  of  coordinates  method  on  the  UTM  projection  surface.  Scale 
factors  were  computed  for  each  line.  The  program  is  given  in  the  Appendix  1  and 
can  be  followed  by  comments  provided  in  the  code. 

The  measured  distances  from  the  system,  for  each  case,  are  corrected  for 
their  slope  with  the  assumption  of  a  right  triangle  (ignoring  the  curv  ature  of  the 


Table  3.  REVISED  POSITIONS  OF  MINI-RANGER 
STATIONS 


Name 

Laiitude(N) 
Longitude(W) 
deg  min  sec 

UTM  coordinates 

Elevation 

Trevor 

36  58  28.785 

122  02  31.536 

585260.158 

4092490.287 

51.1 

Wats 

36  54  26.896 

121  50  39.821 

602947.636 

4085231.429 

22.6 

Packard 

36  49  44.404 

121  46  04.911 

609863.147 

4076611.241 

33.4 

Hays 

36  38  33.824 

121  47  45.894 

607621.298 

4055915.205 

136.7 

Hank 

36  36  25.393 

121  55  08.683 

596670.450 

4051826.796 

134.1 

Earth).  This  assumption  inserts  an  error  of  4  to  5  cm  at  the  extreme  ranges  which 
is  insignificant  for  Mini-ranger  measurements  in  the  areas  covered  by  the  net¬ 
work.  but  it  makes  the  calculations  simpler.  For  the  complete  correction  one  can 
look  in  chapter  15  of  Laurila  [1981]. 

The  next  correction  that  was  made  to  the  measured  distances  was 
for  the  scale  factor  of  the  UTM  projection.  The  following  formula  computes  the 
scale  faetor  for  the  line  between  the  Master  and  each  Reference  station: 

,  El  ^ 

m  =  + - - ) 

6R^ 
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Where  m  is  the  required  scale  factor,  is  the  scale  at  ‘he  central  meridian 
(0.9996),  R  is  the  mean  radius  of  the  Earth  (6375  km)  and  is  the  mean  Easting 
value  for  each  case.  The  mean  Easting  value  is  derived  from  the  formula: 

El  =  e\  +  £,£,  +  £| 

Where  £)  and  £j  are  the  true  Eastings  (x,  —  500000),  of  the  Master  and  the  Ref¬ 
erence  stations  respectively.  The  formula  is  a  simplification  of  the  one  given  by 
Bomford  [1980,  pi 93]  and  the  error  due  to  this  is  of  millimeter  order,  which  for 
our  purposes  is  negligible. 

The  outputs  of  the  program  are  the  least  squares  estimate  for  the  position 
of  the  vessel,  the  variance  covariance  matrix  Z,,  the  residual  \ector  (V)  and  the 
a  posteriori  variance  of  unit  weight  cl.  Measurements  were  assumed  independent 
so  the  diagonal  weight  matrix  is  treated  as  a  vector  in  the  calculations. 

The  a  posteriori  variance  of  unit  weight  cl  is  the  one  that  is  used  for  as¬ 
sessing  both  the  accuracy  of  the  network  and  the  ranging  data,  cl  can  be  scaled 
as  is  well  known,  by  increasing  the  standard  de\  iation  of  a  single  range  meas¬ 
urement.  The  resulting  cl  is  examined  statistically  with  the  chi  squared,  two 
tailed  test,  at  a  95”  o  confidence  level.  The  null  hypothesis,  (?;=  1,  was  tested 
against  the  alternati\e  cr*  #  1.  When  the  null  hypothesis  is  accepted,  it  is  an  in¬ 
dication.  assuming  no  data  blunders,  that  the  system  works  within  the  specifica¬ 
tions  pro\  ided  by  the  manufacturers.  When  the  null  hypothesis  is  rejected  se\  eral 
possible  reasons  for  this  exist.  These  were  summarized  in  section  111- A. 

It  also  should  be  noted,  however,  that  the  presence  of  any  unresohed 
systematic  effects  tends  to  invalidate  this  type  of  test.  One  potential  source,  not 
modelled  in  the  position  computation  algorithm,  involves  the  non-simultaneity  of 
range  measurements  on  a  mo\ing  ship.  In  effect,  the  antenna  is  not  stationary 
over  a  cycle  of  range  measurements  to  all  Reference  stations.  Given  each  range 
measurement  takes  about  55  milliseconds,  four  range  measurements  may  be 
spread  over  a  period  of  about  0.15  seconds.  At  a  nominal  vessel  speed  of  6  knots, 
the  possible  error  is  about  0.45  cm.  which  is  of  little  importance  in  this  analysis. 
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Additional  error  can  be  inserted  due  to  mast  motion.  In  rough  sea  condi¬ 
tions  the  error  is  big  enough  (more  than  1  meter)  and  it  has  to  be  taken  into  ac¬ 
count  by  increasing  the  estinrated  standard  deviation  of  the  range  measurements. 

2.  First  data  set 

The  results  of  the  chi  square  test  at  a  95%  confidence  level  using  standard 
deviations  of  2  and  3  meters  on  the  measured  ranges  respectiv  ely,  appear  in  Table 

4. 

The  acceptable  values  of  ol  under  these  conditions  are: 

For  three  ranges  0.00  <  al  <  5.024 

For  four  ranges  0.025  <  cl  <  3.69 

The  numbers  given  in  the  first  two  columns  of  Table  4  represent  the  per¬ 
centage  of  the  number  of  positions  for  which  the  null  hypothesis  was  rejected 
against  the  total  number  of  the  positions  c.xamined.  The  number  provided  does 
not  include  some  abnormal}'  high  a  posteriori  variances  of  unit  weight  which  were 
caused  by  multipath  effects,  or  low  SS.  These  data  were  rejected  as  blunders. 
Their  percentage  appears  in  the  last  column. 

As  can  be  seen,  even  when  a  =3m,  the  number  of  rejections  is  clearly 
much  higher  than  the  5‘'b  which  one  would  normally  e.xpect. 

In  the  column  'Difference'  we  can  see  how  many  of  the  positions  that 
were  rejected  when  a  =2m.  are  accepted  with  c  =3m.  This  is  an  indication  ei¬ 
ther  of  poor  calibration  of  the  system  or  that  the  system  does  not  perform  to  the 
standard  claimed  by  the  manufacturer.  It  should  be  noted  that  the  range  data 
had  been  corrected  for  an  estimate  of  systematic  errors  determined  in  a  prior 
calibration. 

The  results  shown  in  Table  4  suggest  the  possibility  that  the  master  sta¬ 
tion  used  may  have  beeen  poorly  calibrated.  For  this  reason,  an  examination  of 
the  distribution  of  the  residuals  for  each  reference  station  was  undertaken.  After 
examining  the  distribution  of  the  range  residuals  for  each  reference  station,  it 
became  clear  that  some  stations  were  consistently  biasing  the  position  solutions. 
From  the  distribution  of  these  residuals,  the  bias  was  estimated  by  comparing 
their  mean,  median  and  most  probable  values.  The  mean  residuals  and  the 


Table  4.  REJECTED  PERCENTAGE  WITHOUT  REVISED 
CORRECTORS 


AREA 

B 

B 

Difference 

Rejected 

blunders 

A 

68% 

48% 

20% 

11% 

B 

82% 

37% 

45% 

9% 

E 

69% 

22% 

47% 

25% 

F 

27«'o 

13®'o 

14% 

33% 

G 

43«o 

27% 

16% 

29% 

adopted  correctors  are  shown  in  Table  5.  This  simple  averaging  of  obserx  ed  res¬ 
iduals  is  not  an  ideal  procedure  for  estimating  bias  errors.  The  least-squares 
process  tends  to  redistribute  an  error  in  any  one  range  to  errors  (residuals)  in  all 
ranges.  Recognizing  this  limitation,  however,  the  procedure  gives  a  useful  esti¬ 
mate  of  uncorrected  bias  errors. 

The  results  for  the  corrected  data  are  shown  in  Table  6  and  can  be  di¬ 
rectly  compared  with  Table  4.  The  improxement  as  we  see  is  aboxe  65  ®o. 

In  the  column  'Difference'  we  can  see,  in  comparison  xvith  the  one  of  table 
4.  that  the  x  alues  are  much  lower,  xx'hich  is  an  indication  that  most  of  the  cali¬ 
bration  uncertainty  has  been  remoxed. 

Exen  after  this  procedure,  some  a  posteriori  xariances  of  unit  xvcight  are 
still  high,  especially  in  areas  A  and  B.  The  possible  reasons  for  this  can  be: 

•  Incomplete  remox  al  of  the  blunders. 

•  Operation  of  some  stations  at  loxv  SS  xvhich  in  turn  causes  larger  dex  iations 
from  the  normal  values. 

•  Incomplete  remoxal  of  the  bias  errors  described  aboxe.  The  technique  used 
does  not  eliminate  the  entire  error,  but  only  estimates  it. 

•  For  any  one  of  a  number  of  reasons  the  signal  can  be  lost  temporarily. 
When  it  is  regained  it  sometimes  lacks  stability  for  scxeral  readings. 
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Table  5.  MEAN  RESIDUAL  AND  VALUES  USED  FOR 


RANGE  CORRECTORS 


STATIONS 

Mean  residual 

Used  corrector 

Trevor 

3.91  m 

-3.5  m 

Wats 

-1.66  m 

1.5  m 

Packard 

-0.24  m 

0.0  m 

Hays 

2.61  m 

-2.5  m 

Doppler 

3.64  m 

-3.5  m 

Hank 

-0.24  m 

0.0  m 

Table  6.  REJECTED  PERCENTAGE  AND  IMPROVEMENT  WITH  RE¬ 


VISED  CORRECTORS 


AREA 

0  =  1 

<7  =  3 

Difference 

Improxc- 

ment 

A 

38«^o 

22"  0 

16"o 

44"  0 

B 

31% 

l4"o 

17"'o 

62  "o 

E 

14% 

06"  « 

8"o 

79  "  » 

F 

S-’o 

0"o 

8"o 

70"  ■« 

G 

13% 

7"'o 

6"o 

69"  0 

•  Higher  N  alue  of  the  standard  deviation  of  the  system  used  than  the  one  used 
for  the  weight  matrix  formation. 
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3.  Second  data  set 

The  second  data  set  was  examined  statistically  in  the  same  way  and  the 
results  appear  in  Table  7,  We  can  see  that  the  second  data  set  gh  es  much  better 
solutions.  This  was  expected  because  it  was  known  that  the  master  station  used 
for  this  data  collection,  was  well  calibrated. 

Table  7.  REJECTED  PERCENTAGE  FOR  THE  SECOND 
DATA  SET 


The  pre\  ious  data  set.  for  reasons  noted  earlier,  showed  major  problems 
when  the  ship  was  in  areas  A  and  B.  The  second  data  set.  howc\er.  show^  con¬ 
sistent  results  in  all  areas  and  a  low  rejection  percentage  e\en  when  o  =2m. 
These  results  suggest  that  the  standard  deviation  of  the  range  measurement 
should  be  between  2  and  3  meters. 

E.  GEOMETRICAL  ACCURACIES 

The  position  quality  given  from  the  Mini-ranger  network  in  Monterey  Bay 
depends  strongly  on  the  relati\e  geometry  of  the  ship  and  the  reference  stations. 
The  geometric  position  quality  can  be  represented  in  terms  of  error  ellipses,  con¬ 
fidence  circles,  or  with  several  other  appropriate  techniques.  When  se\eral 
stations  are  used,  computer  techniques  are  necessary  if  detailed  analysis  is  to  be 
undertaken.  Such  techniques  are  described  in  Perugini  [1988]. 
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The  provided  figures  in  Appendix  2  were  taken  from  computer  calculations 
using  software  associated  with  MBARTs  Navigation  processor.  The  numbers 
provided  represent  the  semi-major  axis  of  the  predicted  error  ellipse  at  95%  con¬ 
fidence  level  that  is  formed  at  the  particular  position,  assuming  the  standard  de¬ 
viation  of  a  single  range  measurement  is  3  meters. 

The  best  working  combination  appears  in  Figure  7  and  is  obtained  using 
stations  TREVOR,  PACKARD,  HAYES  and  HANK.  We  can  see  that  the  ge¬ 
ometry  of  the  position  is  very  good  with  two  combinations  of  only  three  stations, 
i.e;  station  combinations  TREVOR,  PACKARD,  HANK  and  TREVOR, 
PACKARD,  HAYES  (see  Appendix  2) 

As  is  easily  seen,  station  TREVOR  occupies  a  position  that  gites  strong  ge¬ 
ometry  in  the  mid- Bay  regions  and  its  position  is  the  most  critical  for  the  network. 
On  the  North  side  of  Monterey  Bay  there  are  also  other  higher  positions  on  which 
a  permanent  station  could  be  installed  with  a  much  better  view  of  the  Bay.  Most 
of  them  are  in  the  Unixersity  of  California  Santa  Cruz  campus  and  some 
arrangments  ha\e  to  be  made. 

In  the  Watsonville  area  the  hills  over  Buena  Vista  Point  ha\e  a  good  Bay 
\iew  but  the  problem  of  a  permanent  electrical  supply  and  security  has  to  be 
soh  ed.  Such  a  position  would  overcome  the  height  problem  of  the  W.A.TS  station 
while  maintaining  its  geometry. 

Stations  PACKARD  and  HAYS  are  at  the  highest  elevations  in  their  areas. 
There  is  no  adx  antage  to  shifting  their  positions  without  deteriorating  the  geom¬ 
etry. 

Station  DOPPLER  does  not  add  much  to  the  geometric  quality  of  the  net¬ 
work  in  the  mid-Bay  but  can  be  used  for  positioning  near  Monterey  harbor  where 
signals  from  PACKARD  and  HANK  arc  obstructed.  It  also  could  be  shifted  NW 
of  Point  Pinos  on  the  lighthouse  or  on  another  high  position  to  improxe  the  ge¬ 
ometry  in  the  Monterey  canyon  area  Further  inxestigation  on  this  station  has 
to  be  done  for  the  best  choice  of  a  permanent  position. 
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Figure  7.  BEST  WORKING  COMBINATION  OF  THE  MINI-RANGER 
STATIONS 
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Station  HANK  gives  very  good  geometry  but  its  position  has  some  limitations 
because  of  obstructions  to  its  Bay  view.  Some  station  at  the  southern  site  of  the 
Bay  is  necessary  to  maintain  a  good  network  geometry  for  the  areas  of  interest. 

Generally  we  can  say  that  the  geometry  of  the  network  is  not  responsible  for 
the  big  standard  deviations  observed  in  most  regions  around  the  Bay,  but  it  be¬ 
comes  more  important  as  one  moves  offshore  near  the  Northen  and  Southern  re¬ 
gions, 

F.  ATMOSPHERIC  CONDITION  EFFECT  OF  THE  AREA  ON  THE 
MICROWAVE  PROPAGATION  VELOCITY 

An  int  estigation  was  carried  out  to  see  the  maximum  possible  error  due  to 
propagation  \elocity  which  could  occur  under  extreme  atmospheric  Bay  condi¬ 
tions  as  compared  with  mean  Bay  conditions.  The  atmospheric  formulation  is 
described  in  Chapter  6  p.l33  of  Laurila  [1981].  The  effects  of  temperature  (T), 
pressure  (P)  and  humidity  (E)  on  the  refractixe  index  (N)  are  gi\en  from  the 
Essen  formula: 


p  I ">  9'>  in 

A-  =  77.62  Y  -  ( -  37.1 9  "Y 

Where  P  and  E  are  in  mb  and  T  is  the  dry  bulb  temperature  in  K°.  But 

^  {RH)e, 

^  100 

Where  RH  is  the  relatixe  humidity  (Vo)  and  Ej  is  the  saturation  vapor  pressure, 
gix  en  in  mb.  at  the  dry  bulb  temperature.  Ej  is  a  logarithmic  function  of  temper¬ 
ature  so  the  effect  on  the  refractixe  index  depends  on  the  temperature  and  the 
relatixe  humidity  as  shoxvn  in  Figure  8,  for  various  relatixe  humidities  and  con¬ 
stant  atmospheric  pressure  at  1013  mb. 

Daily  xveather  logs,  kept  by  the  Meteorology  Department  of  the  NPS.  xxere 
searched  for  the  maxima  and  minima  of  temperature  and  relatixe  humidity  in 
Monterey  during  1988-89.  This  xvas  done  so  as  to  compute  the  ax  eragc  x  ariability 
of  the  refractixe  index  and  the  possible  error  in  distance  that  is  introduced  for  the 
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Figure  8.  VARIATION  OF  THE  REFRACTIVE  INDEX 


actual  atmospheric  conditions.  The  base  for  the  comparison  is  the  refractixe  index 
value  (N  =  320)  that  is  used  from  the  Mini-ranger  for  the  range  calculation. 

In  Monterey  Bay  during  the  spring,  temperature  axerages  about  55  ‘/"and 
humidity  xaries  from  30%  during  the  hotest  days  to  95%  during  the  cold  nights. 
During  working  hours,  a  relatixe  humidit}’  of  55%  with  temperatures  70  F°  aie 
usual.  Extreme  xalues  of  95%  RH  and  55  ‘’/'arc  rare,  but  if  they  are  accompa¬ 
nied  with  high  atmospheric  pressure,  1030  mb,  then  they  gixc  N  =  343.6  which 
in  turn  corresponds  to  a  distance  error  of  23.6  ppm.  In  a  possible  distance  of 
40000  m  the  introduced  error  is  about  1  meter.  On  the  other  hand,  loxv  atmo- 
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spheric  pressure  of  1000  mb  with  80  "F  and  20%  humidity,  gives  N  =  288  and 
32  ppm  error,  which  in  distance  terms  means  1.28  m  at  40000  m  measured  dis¬ 
tance.  Axerage  xalues  of  N  for  the  spring  during  working  hours  are  between  310 
and  330,  which  indicates  that  the  assumed  Refractive  Index  of  320  is  fair. 

During  the  summer  the  increased  temperature  is  accompanied  with  a  drop  in 
the  relative  humidity  so  the  result  tends  to  be  the  same  as  in  spring  with  an  ax- 
erage  N  =  320  to  325.  The  extreme  cases  here  also  rarely  introduce  errors  more 
than  I  meter  at  distances  of  40000  m. 

The  fall  season  is  the  one  with  a  slightly  different  mean  N  and  more  vari¬ 
ability.  Warm  xveather  xvith  high  humidity  is  more  likely  to  occur.  In  the  extreme 
cases,  at  maximum  distance,  the  error  could  reach  2  meters,  but  it  rarely  occurs. 
The  axerage  N  =  330  is  better  to  use. 

Winter  shifts  the  refractixe  Index  back  to  the  loxver  values  due  to  the  loxver 
temperatures  of  the  atmosphere,  xvhich  makes  the  extreme  cases  more  rare.  A 
mean  x  alue  of  N  =  315  is  better  for  use. 

Generally  the  resulting  error  should  be  less  than  a  meter  at  extreme  ranges 
and  sexeral  decimeters  on  axerage,  xvhen  e.xtreme  atmospheric  conditions  happen. 
Exen  though  the  error  is  small,  the  use  of  the  seasonal  refractixe  index  is  adx  ised. 

G.  RANGE  HOLES 

Range  holes  occur  due  to  the  xertical  multipath,  xvhich  is  caused  from  re- 
fiection  of  the  signal  from  the  sea  surface.  This  is  common  to  all  microxvaxe  sys¬ 
tems  that  use  a  line-of-sight  propagation  path.  Figure  9  [Gilb-Weedon.  1976] 
illustrates  the  situation. 

We  count  the  range  holes  from  the  first  to  the  Nth.  Here  N  represents  the 
integer  number  of  xvaxelengths  difference  between  the  direct  and  the  reficctcd 
paths.  The  first  range  hole  occurs  at  the  longer  distance  from  the  shore  station. 

Sexeral  factors  affect  the  range  hole  characteristics: 

•  The  distance  between  the  reference  and  the  master  stations.  The  longer  the 
distance,  the  xvider  the  range  hole  zone. 

•  The  height  of  the  stations.  The  higher  the  reference  stations  the  fewer  the 
range  holes  near  them. 


•  The  wavelength  of  the  signal.  This  affects  both  the  width  of  the  range  holes, 
and  the  distance  from  the  reference  station  at  which  they  occur. 

•  The  atmospheric  conditions  along  the  signal  path  cause  loss  of  SS.  Thus  the 
reflected  signal  from  the  sea  surface,  due  to  additional  losses  from  the  re¬ 
flection.  is  too  weak  to  interfere  the  direct  one,  when  they  both  arri\e  to  the 
master  station. 

•  The  sea  surface  roughness.  This  affects  the  type  of  reflection  which  occurs. 
In  a  strong!}'  wa\e  modulated  sea,  the  appearance  of  the  range  holes  is  rare. 

•  Changes  in  the  sea  le\cl  due  to  tides,  winds  etc.,  cause  range  holes  to  shift 
location  horizontally. 

Small  reflection  angles  cause  more  frequent  range  holes  because  the  signal  has 
lesser  losses  in  comparison  with  the  case  when  the  reflection  angle  is  bigger.  The 
additional  losses  of  SS,  in  the  second  case,  weaken  the  reHecteJ  signal  which  can 
not  then  interfere  with  the  direct  one. 

In  reality,  the  desiructixe  interference  is  rarely  total  due  to  losses  of  the  le- 
flected  signal  from  the  sea  surface,  but  it  can  reduce  the  SS  to  below  the  sensi- 
ti\  ity  of  the  receiver  or  below  the  CST  in  the  case  of  Mini-ranger.  This  will  result 
in  signal  loss,  or  worse,  in  a  bigger  standard  de\ialion  of  the  measurement. 
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In  the  Monterey  Bay  network,  the  range  hole  predictions  are  shown  in  Table 
8  for  each  station.  For  the  maximum  range  of  each  station  we  assume  a  ship  an¬ 
tenna  of  10  m. 

Range  holes  closer  than  the  third  occur  rarely  because  of  the  bigger  angle  of 
incidence  which  attenuates  more  of  the  reflected  signal. 


Table  8.  RANGE  HOLE  PREDICTION  FOR  MONTEREY  BAY 


STATIONS 

Maxi¬ 

mum 

range 

Range  for  the 
first  hole 

Range  for  the 
second  hole 

Range  for  the 
third  hole 

Trevor 

42000  m 

18972  m 

9486  m 

6324  m 

Wats 

30000  m 

8477  m 

4238  m 

2826  m 

Packard 

37000  m 

12440  m 

6220  m 

4147  m 

Hays 

61000  m 

50348  m 

25174  m 

16783  m 

Doppler 

26000  m 

3670  m 

1835  m 

1223  m 

Hank 

60000  m 

1 

49358  m 

1 

24679  m 

16453  m 

Table  8  is  based  on  the  following  formula  [Gilb-Weedon,  1976].  in  which  a 
plane  earth  is  assumed.  R  represents  the  distance  from  the  reference  station  to 
the  range  hole  center  in  meters. 


R  =  2h^ 


!!i 

kl 


The  distance  of  the  ma.ximum  expected  signal  is: 


i2k+\) 


Where  //,  and  lu  arc  the  heights  in  meters  of  the  reference  and  the  master  stations 
respccii\cl\'.  1  is  the  wa\clcngth  and  k  is  an  integer  (k  =  1.2.3....;. 


In  the  Bay  there  is  a  sea  level  variation  due  to  tidal  effects  of  6  to  7  feet 
maximum.  This  effect  tends  to  shift  the  range  hole  center  towards  the  shore  sta¬ 
tion  at  high  water  and  away  from  the  shore  station  at  low  water.  The  o\crall  ef¬ 
fect  of  this  phenomenon  for  the  stations  appear  in  Figure  10.  for  each  range  hole. 

Typical  wind  and  wave  conditions  in  the  Bay  are  such  that  'perfect'  signal 
reflections  are  not  expected.  Thus  range  hole  problems  should  be  rare  and  limited 
to  conditions  of  calm  sea  accompanied  by  light  winds.  Range  hole  phenomena 
were  not  detected  in  the  positioning  data  studied. 
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H.  CRITICAL  SIGNAL  THRESHOLD  (CST)  DERIVATION 

The  relatively  large  distances  and  varying  propagation  conditions  encount¬ 
ered  in  the  Monterey  Bay  network  result  in  many  low  signal  strength  readings. 
It  is  important  that  a  method  be  found  for  using  these  readings  if  at  all  possible. 
Casey,  1982,  has  shown  that  low  signal  strength  ranges  can  be  used  if  corrected 
in  postprocessing  by  an  empirically  derived  relation  between  bias  error,  standard 
deviation,  and  signal  strength. 

The  recommended  procedure  to  establish  this  relation  is  to  perform  repeated 
range  measurements  o\er  a  known  length  baseline,  with  a  variable  attenuator  to 
artificially  \’ary  the  receit’ed  signal  strength.  These  attenuators  were  not 
a\  aliable.  but  the  effect  was  achie\ed  by  ranging  over  relatively  long  baselines  of 
se\eral  different  lengths,  in  order  to  obtain  readings  over  a  spread  of  low  and 
high  signal  strengths.  Two  sets  of  range  readings  were  collected  for  different 
combinations  of  Master  and  icference  stations,  with  signal  strengths  varying  from 
about  5  to  25. 

For  each  Master  Reference  station  pair,  ranges  with  signal  strength  over  15 
were  a\eraged  to  estimate  a  mean  range,  free  of  low  SS  effects.  Differences  be¬ 
tween  this  mean  range  and  all  low  SS  (  <  15)  ranges  were  computed,  representing 
the  additional  error  due  to  low  SS.  By  this  technique,  error  estimates  from  dif¬ 
ferent  Master  Reference  station  pairs  could  be  combined  without  concern  about 
systematic  differences  between  equipment  pairs.  The  resulting  range  error  esti¬ 
mates  were  grouped  by  signal  strength;  5-6,  7-8,  9-10  and  so  forth.  Mean  errors 
and  standard  de\iations  were  then  computed  for  each  signal  strength  group  as 
shown  in  Table  9.  As  expected,  both  the  mean  error  and  standard  deviation  tend 
to  increase  as  the  signal  strength  decreases. 

The  two  data  sets  were  compared  to  ensure  repeatability  of  the  results  before 
deri\  ing  the  empirical  relation.  In  Table  9  we  can  sec  the  results  of  the  statistical 
tests  for  the  two  data  sets.  They  are  treated  as  samples  from  one  population,  but 
each  with  their  own  estimated  means  and  standard  dc\  iaiions  and  examined  with 
the  Fisher  -  Behrens  statistical  test  (Hamilton.  1964]  at  a  95'‘o  confidence  le\cl. 
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as  follows.  The  T  test  limit  was  taken  from  the  Student's  t  distribution  tabic  with 
degrees  of  freedom  calculated  by: 


dof = 


2  2 
'  «i  ^  ' 


ii[Vh  -  0 


»2("2  -  U 


The  calculated  value  for  the  test  is  given  by: 


(■Yi  -  .X2)  -  (a<i  -  ^2) 

o]  nl  ± 


Where  fff  and  a]  are  the  estimated  variances,  yu,  and  are  the  estimated  mean 
errors,  and  «,  and  are  the  number  samples  in  the  first  and  the  second  data  sets 
respectively. 

We  examine  the  hypothesis: 

H^:x\  -  =  0  The  two  data  sets  belong  to  the  same  population. 

Against  the  allernalive 

—  J2  d  The  two  data  sets  don't  belong  to  the  same  population.  Tabic 
9  shows  the  dcri\  ed  statistics  from  the  two  data  sets  and  the  limit  of  the  T  test  in 
each  case.  Every  calculated  value,  for  the  Fisher-Behrens  statistical  test,  that  is 
bigger  than  the  T  \  alue  means  that  the  hypothesis  is  rejected  so  something  else 
going  on  in  the  data  sets. 

In  the  column  'Calculated  value',  the  number  is  derived  using  the  means  and 
the  standard  deviation  of  tin;  measurements.  We  can  see  that  when  SS  is  5-6  and 
17-18,  the  calculated  value  is  bigger  than  the  T  test  limit,  which  means  that  in 
these  cases  the  null  hypothesis  is  rejected. 

The  number  giNen  in  the  column  'Calculated  \alue  with  increased  cr  '  is  a 
combination  of  the  standard  dexiation  of  the  long  term  variation  in  the  instru¬ 
ment  and  its  observed  standard  dexiaiion.  As  we  can  see  the  result  is  a  total 
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Table  9.  STATISTICAL  TESTS  FOR  THE  TWO  SAMPLE  DATA  SETS 


s.s 

1st  data  set 
(mean  error  /  a) 

2nd  data  set 
(mean  error/  a  ) 

T  test 

limit 

Calcu¬ 
lated 
value  (u) 

Calcu¬ 

lated 

value  with 

increased 

a 

5-6 

15.1,4.5 

13.5/5.3 

1.98 

2.48 

1.78 

B 

8. 1/3.7 

7.8;3.8 

1.98 

0.63 

0.41 

9-10 

4.5/3.0 

3.8  2.6 

1.98 

1.68 

0.98 

11-12 

1.9  1.7 

1.9/ 1.7 

1.98 

0.28 

0.19 

13-14 

0.3  1.8 

0.3, 0.7 

1.98 

0.28 

0.2 

15-16 

o 

1 

0.1, 0.6 

1.98 

0.98 

0.37 

17-18 

0.8  1.4 

0.3  0.7 

1.98 

3.35 

1.17 

statistical  agreement  of  the  data  sets.  W’e  conclude,  therefore,  that  they  belong  to 
the  same  population. 

The  next  step  is  the  deri\ation  of  the  equation  that  will  permit  the  use  of  the 
lower  SS  more  erfecti\ely.  The  equation  deri\ed  was 


T 


.  71.093 
^  Inx 


-27.206) 


Where  y  represents  the  mean  error  in  meters  and  x  the  low  signal  strength 
and  the  given  numbers  represent  meters.  The  resulted  y  value  represents  the  ob¬ 
served  minus  the  calculated  (C-O)  range  correction. 

Figure  13  shows  this  resulting  best  fitted  cur\e.  The  correlation  for  the  case 
is  more  than  0.999  which  indicates  that  the  resulting  equation  is  tcry  close  to 
perfectly  modelling  the  data.  This  equation  can  be  included  in  the  software  as  a 
correction  for  the  lower  SS  rather  than  rejecting  them  below  the  CST. 
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Figure  1 1.  FITTED  CURVE  TO  THE  COMBINED  DATA  FOR  TESl  S  I  AND 
2 

Addilionally,  due  to  the  increased  standard  dct  ialion  of  the  ranges  at  lower 
SS  a  higher  value  for  standard  deviation  of  a  single  measurement  should  be  used 
in  the  weight  matrix  formulation. 

I.  VARIATION  OF  THE  STANDARD  DEVIATION  WITH  TIME 

Eigth  sets  of  range  measurements  were  collected  o\er  the  known  dislaiice  be¬ 
tween  station  HAYS  and  geodetic  station  RANGE  7.  The  choosen  distance  was 
short  (2341  m)  to  eliminate  the  effect  on  distance  of  the  variation  of  wcathci 


conditions.  These  8  sets  were  collected  on  8  different  days  and  they  have  length 
of  about  1  hour,  with  3  second  intervals  between  range  measurements.  Table  10 
gi\  es  the  bias  of  these  measurements  with  respect  to  the  known  distance  and  the 
standard  det  iation  of  a  single  measurement.  The  results  from  Table  10  show  that 
the  system  has  a  standard  deviation  which  is  indicative  of  its  precision  rather 
than  its  accuracy.  In  addition,  it  has  a  variable  bias.  The  total  standard  deviation 
for  the  system  must  reflect  both  these  components. 


Table  10.  STANDARD  DEVIATION  AND 
BIAS  OF  RANGE  ERROR 


NUM 

Standard  de\  i- 

ation 

Bias 

m 

0.285 

-0.725 

2 

0.206 

1.447 

3 

0.255 

1.324 

n 

0.266 

1.531 

5 

0.374 

0.577 

6 

0.856 

0.007 

7 

0.576 

0.501 

8 

0.541 

1.173 

The  biases  shown  in  Table  10  were  used  statistically  to  show  that  the  stand¬ 
ard  de\  iation  of  one  measurement  is  2.54  meters,  which  concurs  with  the  \alue 
deri\  ed  by  e\  aluation  of  c?;  estimates  in  section  C  of  this  chapter.  The  procedure 
followed  for  this  deri\  ation  was  a  simple  multiplication  of  the  standard  de\  iation 
of  the  mean  bias  found  from  the  data  set  by  the  square  root  of  the  number  of 
data  sets.  The  user  shouldbe  aware  that  the  standard  deviation  can  be  larger  on 
board  a  \essel  in  a  rough  sea  due  to  mast  mo\emenls  and  the  selecti\  iiy  of  the 
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stations  that  the  system  makes(non  simultaneous  measurements  to  the  reference 
stations).  The  bias  variation  from  day  to  day  makes  the  calibration  of  the 
stations  time  consuming  because  a  one  day  calibration  will  not  remot  e  all  the 
bias. 


Table  11.  MEAN  DIFFERENCES  AND  STANDARD  DEVIATIONS  FOR 
STATIONS  HAYS  AND  TREVOR 


DIS¬ 

TANCE 

Mean  differ¬ 
ence  (Hays) 

Standard  devi¬ 
ation  (Hays) 

Mean  differ¬ 
ence  (Trexor) 

Standard  de¬ 
viation 
(Trexor) 

SHORT 

0.55 

0.8 

- 

- 

ME¬ 

DIUM 

1.02 

0.55 

0.91 

0.48 

LONG 

1.27 

0.75 

1.48 

0.68 

V.LONG 

1.28 

1.2 

0.84 

0.95 

J.  VARIATION  OF  THE  STANDARD  DEVIATION  WITH  DISTANCE 

The  second  thing  was  the  variation  of  the  standard  deviation  with  the  dis¬ 
tance  or  with  S.S,  when  it  is  above  the  CST.  Several  problems  limited  the  data 
collection,  like  the  unavailability  of  an  attenuator,  the  poor  network  of  known 
positions  around  the  Bay  and  the  low  elevation  of  possible  calibration  sites.  The 
data  were  sufficient  only  for  the  evaluation  of  station  HAYS  and  partially  of 
station  TREVOR. 

Table  11  shows  the  mean  difference  of  the  measured  distances  from  the  true 
distances  and  the  standard  dex  iations  of  each  set  of  range  measurements.  We  can 
see  that  there  is  variation  of  the  mean  error  with  distance  in  both  stations  but  no 
indication  exists  that  a  function  can  be  dcrixed  because  the  mean  differences  with 
the  gixen  standard  dex  iations  are  not  statistically  different.  Actually,  a  functional 
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variation  of  the  standard  det  iation  with  distance  was  not  expected,  when  oper¬ 
ating  above  CST. 
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V.  DISCUSSION  AND  CONCLUSIONS 

The  Mini-ranger  network  in  the  Monterey  Bay  was  installed  for  research 
purposes  by  MBARI  and  NFS  and  was  examined  from  various  aspects.  Useful 
results  were  drawn  that  can  help  in  the  understanding  of  the  operation  of  trans¬ 
ponder  type  systems  in  other  areas,  and  that  will  make  the  present  network  more 
useful  when  operating  at  low  SS. 

Figure  12  shows  the  overall  coverage  of  the  Mini-ranger  with  all  the  limita¬ 
tions  taken  into  acount,  except  that  weather  conditions  can  attenuate  the  signal 
and  in  turn  reduce  the  maximum  range  of  the  stations.  For  the  maximum  range 
derivation  a  height  of  the  master  station  on  the  ship  was  assumed  to  be  10  meters. 

From  this  network  and  the  data  processed,  the  following  conclusions  can  be 
drawn: 

•  For  the  calibration  of  such  systems,  if  one  wants  the  ma.ximum  accuracy  that 
can  be  proA  idcd,  a  bias  taken  from  single  obser\'ations  over  known  distances 
or  with  the  baseline  method  may  be  not  enough.  It  is  better  to  use  repeated 
observations  and  a  mean  corrector  derivation. 

•  The  positions  of  the  stations  in  the  Bay  were,  with  the  e.xception  of  station 
HANK,  determined  to  a  sufficint  le\  el  of  accuracy  such  that  they  were  not 
responsible  for  some  of  the  large  range  standard  de\  iations. 

•  The  software  of  the  user  is  important  in  the  position  derh  ation  and  can 
cause  unexpected  errors  if  not  properly  validated. 

•  The  meteorological  conditions  in  the  Bay  can  cause  attenuation  of  the  signal 
and.  in  turn,  reduce  the  maximum  range  ability.  The  range  error  from  the 
delay  of  the  signal  due  to  meteorological  conditions  does  not  exceed  the  1.2 
meters  in  the  extreme  conditions  and  is  smaller  than  0.2  meters  under  the 
usual  operational  conditions. 

•  Range  holes  in  the  Bay  can  occur,  but  should  be  rare  due  to  usual  windy 
Bay  conditions  and  the  modulated  surface  from  the  swell. 

•  We  ha\e  shown  that  ranges  measured  at  low  SS  need  not  be  rejected,  but 
can  be  corrected  using  an  appropriate  error  function.  This  function  may  \  ary 
from  system  to  system. 

•  We  show  that  the  standard  de\iation  can  be  considered  to  consist  of  two 
components,  the  unpredicted  bias,  which  is  smaller  when  the  system  is  well 
calibrated  and  the  measured  standard  de\iation  of  the  data.  Generally  the 
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IN  MONTEREY  BAY 


ideal  conditions  is  a  little  greater  than  2  meters. 
3  meters  when  the  system  operates  at  sea. 
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•  There  are  no  indications  for  variation  of  the  standard  deviation  with  dis¬ 
tance  (or  SS)  above  CTL. 
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APPENDIX  A.  COMPUTER  PROGRAM  LISTING 


REM  PROGRAM  "LESS"  FILE:  NICK3.  BAS 

NICK  KRIONERITIS  Date:  8-8-89 

Modified  10-17-89  by  KJS 

************************************************************* 
REM  This  program  calculates  the  optimum  position  of 
a  ship  using  distances  from  known  stations.  It 
serves  the  solution  for  a  miniranger  network  with 
accuracy  better  than  10  centimeter. 

REM  The  subrutine  SOLVERANGX  was  changed  from  CDR  K.  SCHNEBELE 
program  "HYDROPLOT". 


$STACK  9216 

SINCLUDE  "SOLRAN’GX.  SUB" 

SINCLUDE  "NICKMR. SUB" 

ON  ERROR  GOTO  ERRORHANDLER 

Print  "Enter  the  number  of  stations" 

Input  "N?b=",N°o 

Lprint  "Number  of  stations="  ,N?i 
'  DECLARATION  OF  THE  VARIABLES 

'  AP>'a(1):  Approximate  position  of  theship  in  meters  (UTM) 

'  SP^/(2);  Station  positions  in  2D  array  in  meters  (UT.M) 

'  R^Kl):  Ranges  observed  on  the  receiver  in  meters 
'  L.\»(l):  Stations  Eastings  for  the  scale  factor  affection 

'  Xj'KI):  False  Eastings  of  the  stations 
'  y?;(l):  Northings  of  the  stations 

'  Z'-(l);  Elevation  s  of  the  stations  from  MSL 

’  SF-‘(1):  Scale  factor  calculated  for  each  case 

'  W!(l):  Diagonal  Weight  matrix  assuming  indepedence 

'  among  the  measured  ranges 

'  C0DE°o(l):  Index  numbers  of  the  4  stations  in  each  record 
'  AilSta?r'(2):  UTM  Coords  (x,y,z)  of  6  stations  possible  to  use 

'  Vectors  for  S0LVERANG2  but  not  used  in  this  version 

'  CmOs!(l):  Bias  corrections  for  each  of  6  stations 
'  Wgts!(l):  Weight  factor  for  each  of  6  stations 


DIM  AP4‘(  1:  3) 

DIM  DYN-4MIC  SP?K  1;  N%,  1:  3)  ,  R#(  1:  N%)  ,  LXiK  1:  N%) 

DIM  DYNAMIC  X.^i(l:N‘’c)  ,  Y/>(1:N’»)  ,  Z^C  1:  N?o) 

DIM  DYNAMIC  M.-'>(l:N°c)  ,  SF^K1:N%)  .  BRes#(  1:  N?0  ,  WKliN'o) 

DIM  code°o(l:  4)  ,AllSta^/(  1:  6, 1:  3),Cm0s!  (1:  6),Wgts!  (1:  6) 

? 

PRINT  "FREE  ARRAY  &  STACK  AT  STARTUP",  FRE(-l),  FRE(-2) 
CALL  GetStationsCAllSta.'K  )  ,Cm0s!  ( )  ,Wgts!  (  ) ) 

Lprint  CmOs!  ( 1); CmOs!  (2); CmOs! (3); CmOs! (4); "  CmOs" 

CALL  SetFiles 

Print  "Enter  assumed  position" 

Input  "AX="  ,  AP?';(  1 ) 
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Input  "AY=",AP#(2) 

Input  "AZ=",AP#(3) 

LPRINT  "Assumed  Position"; 

For  i%=l  to  3 

Print  "Assumed  position  in  UTM  (X,Y,Z)="  AP#(i%) 
LPRINT  USING  AP//(i%); 

Next  i% 

'  **************'**********ilr*********'/r**ilr***iHWr* 


REM  Aproximations  for  calling  the  subrutine  S0LVERANG2* 
'  ********************************************** 


SFAC!=1.  0 
SHGT^',f=AP^'/(3) 
AP?'K3)=0.d0 
LM#=1.  1 
Rm)'/=6368459 


'Approximations  for  the  subrutine  for  the 
'  corrections  that  have  already  done 

'Convergence  Limit  test  (meters) 

'Mean  radius  of  Earth  in  meters  *  0.9996 


'  ********************Vr************ 


REM: Declaration  of  Weight  matrix  * 

'  *******************^VVr************* 

For  i%=l  to  N% 

W!  (i”^)=.  25 


'Diagonal  matrix  assuming 
' independent  measurements 
Next  i?o 

'  **Vf***VfV,'V-**********Vc*************** 

Beep  2:  CLS 

PRINT  "Put  Printer  online  in  COMPRESSED  PRINT  Mode!" 

PRINT  "  Hit  any  key  when  ready  to  continue” 

WHILE  NOT  INSTAT 
LOCATE  2,37 
PRINT 
WEND 

Print  INKEY$ 

f 

'  LOOP  TO  READ  AND  PROCESS  ALL  DATA  UNTIL  END  OF  FILE  HIT 


'WHILE  NOT  EOF(l) 
For  k“o=l  to  500 


CALL  READNICKMR(sec&,code%(  )  ,R#(  )  ,X#()  ,  Y/K  ) ,  Z#( ) ,  A1  lSta,v( ) ) 

*  **************Vr***Vr**************** 

REM  Correction  for  slope  distance  * 

'  *******Vc*********')V************Vr**** 


For  i%=l  to  N?o 

If  (R#(n)  >  l.dO),  Then 
DHGT#=Z#(i%)-SHGT# 

HGT//=R^i(i%)  2-DHGT/^  2 

IF  HGT#<  O.dO,  Then 

PRINT  "Warning  -  at  time  secS; "  DHGT  exceeds  Range" 
PRINT  "  Hit  any  key  to  continue" 

While  Not  INSTAT 
Delav  0.  2 
Print 
Wend 
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Print  INKEY$ 

Else 

R#(i%)=SQR(HGTy/) 

End  If 

Else 

RiK  i7o)  =0.dO 
End  If 

Print  "Ranges  corrected  for  slope="  R^/(  i%) 
Next  i7o 

'  *****************"***************■*•***** 

REM  Calculation  of  the  scale  factors  * 

AXP,-!4=AP,ii(l) -500000 
For  i%=l  to  N% 


LX;'K  i7o)=XiK  i%)  -500000 

2  +  AXPy^*LX#(  i%)  +  LXT-'Ki^) 
SFy/(i?^)=0.  9996  +  .9996*M#(i%)  /(6*Rm«) 

Print  "Scale  Factors  for  UTM="  SF#(  i?i) 

'  Lprint  "Scale  factors  for  UTM="SF#(i%) 

REM  Correction  for  the  scale  factors  * 

*  ^VVfVfVrVrVrVfsV^V-V*Vr':VVr“V***Vr*iVVf*'cVc*Vf^V*VV?'5V9V?r'^r*:frVf^VVr 

R?'/(i%)=R,-!Ki%)*FFyKi%) 

Print  "Corrected  ranges="  R#(i%) 

'  Lprint  "Corrected  ranges="  R?K  i/o) 

Next  i% 


Vf****-. 


REM  Formation  of  station  position  matrix 


For  i?o=l  to  N°o 
SP/;(i%,l)=X^Ki%) 
P?'^(i'’,,2)=Yf:i(n) 
SP?Ki'o,3)=0.  do 
Next  i°o 


2 


CALL  SolveRangX(AP#(),SP,-'K),R?'K),CmOs!  (),N?e,SFAC!  ,  W!  ( ) ,  LM^r‘  ,BRes, 
sxx!  ,  syy!  ,  sxy! ,  variance! ) 

CALL  RiteSol(sec&,AP/‘( )  ,code?o( )  ,R7‘'/(  )  ,BRes//( )  ,sxx!  ,syy!  ,sxy!  , variance!  ) 

I 

LPRINT  USING  "  i'/";  sec&  AP^KD  AP.::i(2); 

For  i°o  =  1  to  4 

LPRINT  USING  "  UU.i^  +##.#";  code%(i%)  R^Ki°o)  BResfKi°o); 

Next  i°o 

LPRINT  USING  "  +#M.  #  sxx!  syy!  sxy!  variance! 

'  WEND 

t 


Next  k% 

CLOSE  in 
CLOSE  72 
END 

ERRORHANDLER: 

PRINT  "FREE  AT  ERROR" ,FRE( -1) ,  FRE(-2) 

PRINT  "AN  ERROR  TYPE",  ERR,  "HAS  HAPPEN  AT",  ERADR 
END 
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SINCLUDE  "openfile.  sub" 

' *************  GetStations  Subroutine  ************  kjsiOct  89  * 

'  file: NICKMR.  SUB 

SUB  GetStations( Allstate 2), CmOs!  (D.Wgts!  (1)) 

LOCAL  i7o,  s$,  code% 

CLS:  PRINT  "insert  Disk  with  File  ST.  DAT  in  Drive  B: " 

Print  "  Hit  any  key  to  continue  when  ready"; 

While  Not  INSTAT 
Locate  2,37 
Print  "!" 

Wend 

Print  INKEY$ 

CALL  0PENFILE("I",1,"B:  ST.DAT",128) 

For  i%  =  1  to  6 

INPUT  #1,  code%, Allstate i%,l),AllSta#(i%, 2), AllSta#(i%, 3), CmOs!  (i%),_ 
Wgts!  ( i%) ,s$ 

Next  i% 

CLOSE  in 
END  SUB 

************  SetFiles  Subroutine  *******************  kjs:  Oct  89  ******* 

’  file;  NICKMR. SUB 

1 

'opens  input  data  file  of  Miniranger  values(as  #1),  and  opens  output  file 
'  for  the  computed  results  (as  #2) 

SUB  SETFILES 

LOCAL  s$ ,txt$ ,FiName$ 

CLS:  sS  =  "  (Drive:  path  name. ext)  " 

PRINT  "Enter  MiniRanger  data  file  for  input"  +  s$; 

I.NPUT  FiName$ 

CALL  0PENFILE( "l" , 1 ,FiName$ , 128) 

PRINT  "Enter  solutions  file  for  output"  +  s$; 

I.NPUT  txt$ 

CALL  0PENFILE("0",2,txt$,128) 

PRINT  #2,  "Solutions  from  file  ";  FiName$ 

END  SUB 

' Vf**********  ReadNickMR  Subroutine  ******************  kjs;  Oct  89  ****** 
'  file:  NICKMR. SUB 

'specific  to  read  mini-ranger  files  (from  #1)  from  Krionertis  thesis  work 
'with  time  in  secs  &  4  station  numbers  &  4  ranges. 

SUB  READNICKMR(  sec& ,  code?i(  1 )  ,  rang#(  1 )  ,X#(  1 )  ,  Y#(  1 )  ,  ZiK  1 )  ,  A1  lSta(> ( 2 ) ) 

LOCAL  s  $  ,  i% ,  nucode?o( ) 

DIM  DYNAMIC  nucode^o(  1:  4) 

I 

LINE  INPUT  in,  s$ 
sec&=VAL(MID$(s§,8,7)) 

For  i%  =  1  to  4 

nu cod e'’o(  i'o )  =VAL(  M I D$  ( s  $ ,  1 1+4* i% ,  4 )  ) 
rang.^Ki%)=VAL(MID$(s$,22+9*i%,9)) 
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'check  nucode%(),  if  not  same  as  code%()  then  get  new  station  coordinates 
I  for  the  X//(),Yi'K),Z//()  arrays 


If  (code%(i%)  <>  nucode%( i%) ) ,  Then 
code%(i%)  =  nucode%(i%) 

X^Ki%)  =  AllSta?'/(code%(i%),l) 
Y#(i%)  =  AllSta,-'/(code?4(i%),2) 
Z#(i%)  =  Allstate code%Ci%), 3) 

End  If 
Next  i% 


END  SUB 

'***********RiteSol  Subrutine  **************************  kjs:  Oct  89***** 
’  File:  NICKMR.  SUB 


Writes  Miniranger  Solution  and  error  statistics  to  File 


SUB  RiteSol(sec&,AP#( 1) >code%( 1) ,R#(1) ,BRes/K  1) jsxx!  ,syy!  ,sxy!  ,vari! ) 
LOCAL  i% 


PRINT  n,  USING  secSc  AP^KD  AP,-'K2); 

f$  =  "n  MUihhif  +iHhir 
For  i%  =  1  to  4 

PRINT  n,  USING  f$;  code%(i%)  R/Ki%)  BRes/Ki%); 

Next  i?4 

PRINT  n,  USING  "  m.U  HihU  mHHhU"-,  sxx!  syy!  sxy!  vari! 

END  SUB 


'  *Vf*** 


File  =  SOLRANGX. SUB 

*  *  *  *  *  *  SOLVE  RANGE  POSITION  (ver  X)  *  *  *  SUBROUTINE  *  *  * 


Returns  grid  x,y  in  first  2  elements  of  apsni'K ) 


INPUTS: 

apsnj'rO)  =  approximate  position  as  grid  x,y,z 

spsns#  =  station  positions  in  2D  array  as  grid  x,y,z 

(nsta?i,3)  (array  has  nsta%  rows  by  3  columns  for  x,y,z) 
ORS(Knsta?o)=^  observed  ranges  to  each  station  in  spsnsi^  row  order 
rangcor!  (nsta'o)  =  C-0  corrections  for  each  range  reading 
nsta°i  =  number  of  stations/ranges 

scalefac!  =  grid  scale  factor  (grid: true  ratio) 
wgts! (nsta%)  =  weights  for  each  rang^  (1/sigma 
dxlim?r  =  convergence  limit  to  stop  iterating  position 
solution  (given  in  grid  units  -  e. g.  meters) 

OUTPUTS: 

apsn(K3)  =  x,y  position  returned  in  1st  two  elements 
(elevation  of  antenna,  z,  unchanged) 

BRes(Knsta?i)  =  residuals  on  observed  ranges  (ob-computed) 

sxx!  =  variance  in  x  (from  Inv(ATWA)) 

syy!  =  variance  in  y 

sxy!  =  covariance  of  x  &  y 

variance!  =  est  of  sigma-o  (2F’lW)gt ]B/ihi(Sbs-2) ) ) 


$ INCLUDE  "LUDCMP.  BAS" 

S INCLUDE  "LUBKSB.BAS" 

SUB  SolveRangX(  apsn()(  1)  ,spsns#(2)  ,ORS#(  1) , rangcor!  ( 1)  ,nsta°o, scalefac!  , 
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wgts!  ( 1)  ,dxliin#,BRes#(  1)  ,sxx!  ,syy!  ,sxy!  .variance!  ) 
LOCAL  m  )  ,spsn#( )  ,ORa/^Ax//,Ay#,Bi#,ATtf/K )  .ATWBj'/C )  .ATWAj'K  ) ,  Indx%( ) , 
sum#,iteration%,i%,  j%,k%,YCol^K  )  ,ATWAInv#( )  ,xdx%( ) ,  Btemp.-'K ) 
DIM  spsn^Kl:3),  ATWB#(1;2),  ATWA#(  1:  2 , 1:  2)  ,  Indx%(  1:  2)  ,xdx%(  1:  4) 

^  Btempi'/C  1:  4) 

i%=l 

For  j%  =  1  to  nsta?i 

If  ORs#(j%)  >  l.dO,  then 

xdx?i(i%)  =  j%  'xref  to  non-zero  ranges  (range  >  Im) 

INCR  i% 

End  If 
Next 

nranges%  =  i%-l 

DIM  A#(l:nranges%,l:2),  ATW#(  1;  2, 1:  nranges%) 

SolveAgain: 

For  i%  =  1  to  nranges% 

ORa#=ORS^Kxdx%(  i%)  )  +  rangcor!  (xdx?b(  i%) ) 

For  j%  =  1  to  3 

spsn#(j%)  =  spsns#(xdx%( i%) , j%) 

Next  j% 

CALL  RangeLOPCapsn/K  )  jspsn#( )  jScalefac!  , ORa# ,  Ax# ,  Ay^'/ ,  B i# ) 
A#(i%,l)=Ax# 

A#(i%,2)=Ay# 

BRes#(  i?4)=Bi#  'Residual  (o-c)  on  each  range 

Next  i% 

I 

I  Compute  ATW  C2 ,nranges®i| ,  &  ATWB  C2| 

For  j%  =  1  to  2 
sum#=0.  do 

For  i%  =  1  to  nranges% 

ATV,’#(  j%,i'’o)  =  A#(  i®o,  j%)*wgts!  (xdx%(i®;) ) 
sum#  =  ATW#(  j%,  i%)*BReS7#(  i%)+sum# 

Next  i% 

AWBi^'C  j°o)=sum# 

Next  j?o 

'  Compute  ATWA  C2,2| 

For  =  1  to  2 
For  i^o  =  1  to  2 
sum#  =  0. do 

For  =  1  to  nranges°o 

sum#  =  ATW#  ( j  % ,  k'o  )  *  A#  (  k?i ,  i?i  J + s  urn# 

Next  k% 

ATWA#(  j?o ,  i?o)=sum# 

Next  i?i,j?i 

t 

CALL  LUDCMPC  Alk'A#( ) ,  2 ,  Indx%(  )  ,  i%) 

CALL  LL'BKSB(ATWA#(  )  ,2 ,  Indx?.(  )  ,AWB#(  ) ) 

sum#=0. do 

For  i%  =1  to  2 

apsn#(  i”o)=apsn#(  i?o)+ATWB#(  i®o)  'shifts  position  by  dx,dy  which  has 
sum#  =  ATWBi2(  i-^o^um#  'been  returned  as  ATVv'BC  )  from  LUBKSB 


Next  i% 

sunn';*  =  SQR(sum#) 

If  sum#  >  dxlim#.  Then 
INCR  iteration^ 

If  iteration^  >  7,  Then 
Beep  3 

Print  "Position  did  not  converge  --  NO  SOLUTION! 

Print  Using  "Delta  Position  =  +#.  ####meters";  sum# 
sxx!  =  0.0:  syy!  =  0.0:  sxy!  =  0.0:  var!  =  LOG10(sum#) 

EXIT  SUB 
End  If 

GoTo  SolveAgain 
End  If 

1 

'  Compute  variance-covariance  values  from  inverse  of  ATWA 
'  sxx!  =  sigma-x  2 

'  syy!  =  sigma-y  2 

'  sxy!  =  sigma-xy 

ERASE  ATW#,  spsn#  'make  some  room  in  memory 

DIM  YCol#(l:2),  ATWAInv#( 1: 2, 1: 2) 

For  j/b  =  1  to  2 

For  i?i  =  1  to  2  'identity  vector,  YCol,  to  find 

YCol,#(  i?i)  =  0.  do  'inverse  using  back-substitution 

Next  i°o 
YColfK  j°o)=l.  do 

CALL  LUBKSB(ATWA#() ,2,Indx%(),YCol#()) 

For  i%  -1  to  2 

ATWAInv#( i%, j%)=YCol#( i%) 

Next  i°o 
Next  j°o 

sxx!  =  CSNG(ATWAInv#(l,l)) 
syy!  =  CSNG(A-nv'AInv#(2,2)) 
sxy!  =  CSNG(AT\s’AInv#(l,2)) 

'  Compute  estimate  of  unit  vCtftWgtK;)“/nasaiS41^fttts 
'  Note  that  range  residuals  are  in  BRes#( ) 

If  nranges?o  =  2,  Then 

variance!  =  9999  'signifies  no  variance  computable 

Else 

sum#  =0. do 

For  i?o  =  1  to  nranges'o 

sum#  =  Bres#(i%)  2*wgts! (xdx%( i%) )  +  sum# 

Next  i% 

variance!  =  CSNG(sum#/(nranges%-2) ) 

End  If 

t 

If  nranges?o  <  nsta^o,  Then  'Set  residuals  into  proper  channel 

For  i°o  =  1  to  nsta’o 
Btemp#(i?o)  =  0.  dO 
Next  i°o 

For  j%  =  1  to  nranges?4 

Btemp#(xdx%(  j?,,) )  =  Bres#(j%) 

Next 

For  i%  =  1  to  nsta“o 
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If  nranges%  =  2,  Then  Bres#(i%)  =  0.  dO 
Bres#( i%)  =  Btemp#( 1%) 

Next  i% 

End  If 

t 

END  SUB 

I 

*******  RANGELOP  *  *  SUBROUTINE  ********** 

I 

'Computes  A  matrix  components  Ax  &  Ay  and  grid  range 
' for  input  to  multiple-LOP  solution  algorithm. 

'Inputs:  apsn#()  =  approx  position  x,y,z  coordinates  (grid) 
'  spsn#( )  =  station  x,y,z  coordinates 

'  scalefac!=  grid  scale  factor  (point  or  line) 

'  ORa#  =  observed  slant  range 

'Output:  Ax^/  &  such  that  (Ax)dx+( Ay)dy  =  OR-CR 
'  OmCR#  =  Observed  minus  Computed  range  at  apsn// 

'Note:  Observed  range  (OR)  has  been  corrected  for  elevation 
'  differences  and  grid  scale  before  comparing  with  CR. 


SUB  RangeL0P( apsn#( 1) »spsn#( 1) jscalefac! 

0Ra)¥ ,  Ax// ,  Ay# ,  OmCR#  ) 

LOCAL  n.CR/Z.dx/ZO 
DIM  dx#(l:  3) 


For  i^i  =  1  to  2 


dx#(i%)  =  apsn#(i%)-spsn#(i%) 

Next  i% 

CR#  =  SQR(dx#(l)  2+dx#(2)  2) 

Ax#  =  dx#(l)/CR# 

Ay#  =  dx#(2)/CR# 
dx#(3)=0. dO 

OmCRi#  =  0Ra,#*scalefac!  -  CR// 

END  SUB 


file:  OPENFILE.  SUB 


Opens  a  disk  file  using  TBasic  commands.  Preforms  error  checking  and 
allows  recovery  from  common  mistakes. 


INPUTS:  modes$  =  One  character  designation  for  filemode 

'  "O"  =  for  sequential  output  (to  be  written  to) 

'  "l"  =  for  sequential  input  (to  be  read  from) 

'  Csee  Turbo  Basic  Ref,  pg  282,  for  other  choices! 

'  FileNum%  =  file  number  to  be  used  in  read/write  calls 

'  FileName$  =  Drive: /path/name. ext  of  file 

'  RecLength?o  =  record  size  in  bytes,  1  to  32767,  typical  =  128 

SUB  OpenFile(mode$  ,FileNum?i,FileName$  ,RecLength/o) 

LOCAL  ermum^o,  ans$ 
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OpenSequence; 

ON  ERROR  GoTo  ErrorHandling 

OPEN  mode$ ,FileNum%,FileName$ ,RecLength% 

ON  ERROR  GoTo  0 

EXIT  SUB 


ErrorHandling: 


Beep  2 

errnum%  =  ERR 
Select  Case  errnum% 

Case  53,75,76 

Print  "FILE  or  DRIVE  PATH  Not  Found  for  .  .  .  ",  FileName$ 


Input  "Re-Enter  Drive:  path  name,  ext 
Case  68,71 

Print  "DISK  DRIVE  NOT  READY  for  ",  FileName$ 

Case  Else 

Print  "Opening  File  Gave  Error  Number  ",  errnum°i 
End  Select 

Input  "(R)etry  or  (H)alt  Program  ans$ 

If  Ucase$(ans$)  =  "H",  Then 
Resume  HaltSequence 
Else 

Resume  OpenSequence 
End  If 


FileName$ 


HaltSequence: 


'terminates  main  program 


Print  "Program  Halted!  !  !  " 

Print  "Subroutine  OpenFile  with  Error  Number  ",errnum®o 


Print  encountered  at  program  address  ,  ERADR 
END 

END  SUB 

' LU  Back  Substitution 


'  Ref:  Numerical  Recipies,  by  WH  Press,  page  37 


SUB  LUBKSBC  A?'/ ( 2 ) ,  N% ,  INDX?b(  1 )  ,  Bt'K  1 ) ) 
LOCAL  ll?i,i%, j?b,sum,^/ 


in=o 

For  i%=l  to  N% 
ll%=lndx®4(  i°o) 
sum,''/=B?:f(  ll?i) 

B#(ll%)=B#(i?i) 

If  ii°o<>0,  then 
For  j/o=ii%  to  i%-l 

s  um/|i=s  umi^^  -  A#  (  i% ,  j  %  )  *B/f  (  j  %  ) 
Next  j% 

Elself  sum#<>0.  do  Then 
ii%=i% 

End  If 
B#(  i”o)=sum,''/ 


Next  i% 


For  i%=N%  to  1  Step  -1 
sum#=B#(i%) 

If  i%<N%,  Then 

For  j%=iro+l  to  Nro 

sum,^/=sum#-A#(  i%,  j%) 

Next  j% 

End  If 

B#( i%)=sum#/A#( i%, i%) 

Next  i« 

END  SUB 


t 

I 


|LU  Decomposition  Subroutine 


'  Ref:  Numerical  Receipes,  by  WH  Press,  pages  34-36 


SUB  LUDCMPC  kU2 )  ,N°i ,  Indx?i(  1)  ,D%) 

LOCAL  vv#()  ,Dum,-'/,i%,  j%,k%,imax%,aamax?’/,sum,‘'; 
STATIC  nmax?i,  tiny# 


nmax?b=10 
tiny,#=l.  OD-20 

DIM  vv#(  1:  nmax?i) 


For  i%=l  to  N?o 
aamax#=0 
For  j%=l  to  N% 

If  ABS(A#(i%, j%))>aamax#,  Then  aamax#=ABS(A#(i%, j?b) ) 

Next  2% 

If  aamax#=0,  Then 

Print  "Matrix  has  zero  row  #  "  i% 

Exit  Sub 
End  If 

vv#(i%)=  1/aamax#  'save  row  scaling  factor 

Next  i/o 

1 

For  j%=l  to  N% 

If  j%>l,  Then 
For  i%=l  to  j?o-l 
Sum#=A#(  i?o,  j%) 

If  i?4>l  Then 
For  k?i=l  to  i?o-l 
S  um#=S  urn# - A# ( i% , k% ) *A# ( k% , j  % ) 

Next  k% 

A#(  i%,  j%)=Sum,# 

End  If 
Next  i®i 
End  If 

'  continuing  j  loop,  searching  for  largest  pivot 

aamax#=0 

For  i?i=  j%  to  N% 

Sum#=A#C  i/o,  j?i>) 
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If  j%>l,  Then 
For  k%=l  to  j%-l 
Sum#=Sum# -m  i% ,  k% ,  j%) 

Next  k% 

A^Ki%,j%)=Suin# 

End  If 

Dum#=vv#(i%)*ABS(Sum#)  'figure  of  merit  for  pivot 
If  Dum#>aamax#,  Then 
imax%=i% 
aamax//=Dumj'/ 

End  If 
Next  i% 

f 

If  j%<>imax%j  Then 
For  k%=l  to  N% 

Dum#=A#( imax% , k%) 

AiK  imax% , k%)=A#( ,k%) 

A#(j%,k%)  =  Dum# 

Next  k% 

vv#(  imax?4)=vvj'/(  j?i) 

End  If 

1 

Indx?i( j?o)=iniax%  'divide  by  pivot 

If  (j%<>N%),  Then 

If  A/Kj‘’o,j%)=  0,  Then  A#( j%, j%)=tiny# 

Dum#=  1  /  A//  ( j  ,  j  % ) 

For  i%=j%+l  to  N% 

A^K  i% ,  j  % )  =  Af'K  i% .  j  % )  *'DumU 
Next  i% 

End  If 
Next 

t 


' interchange  rows  if  reqd 


'change  parity  of  D% 


If  A,'<(N%,N%)=0,  Then  A#(N%,Nr,)=tiny# 

END  SUB 

*  i^:V*‘>V?V*.V-/fVfVf*VrVfVfVfVfVfyfVfVfyrVfVfVfVf^V'3VVf^V;iV^VVr^V*^V>V*yrVrVf>V'>VVry(r>V7VVrVf^V^V9V'}VVfVfV? 

*  yfVfyf*yfyr*yf'5Vyf*y.'VfVr*Vr*V«"5V’}V**Vr’5Vyryf^VVcVf^VVrV<”}VVfVfyryriVyr'>VVfycyfVryrVrVc->VVrVrVr’>VyfVc 


>v  y?  Vr  V?  yp  y?  y?  * 
yr  y?  yp  yp  >>  y?  yp  yp 


1, 

585260. 161, 

4092490. 284, 

51.7, 

0.  0, 

0. 

25, 

TREVOR 

2, 

602947. 672, 

4085231. 399, 

23.  1, 

0.  0, 

0. 

25, 

WATS 

3, 

609863. 128, 

4076611.  345, 

33.  9, 

0.  0, 

0. 

25, 

PACKARD 

4, 

607621. 289, 

4055915.  264, 

137. 2, 

0.  0, 

0. 

25, 

HAYES 

5. 

600434. 009, 

4051260.  014, 

10.  0, 

0.  0, 

0. 

25, 

DOPPLER 

6, 

596669. 490, 

4051826.  411, 

134. 5, 

0.  0, 

0. 

25, 

HANK 

APPENDIX  B.  GEOMETRICAL  ACCURACIES  OF  THE  IMONTEREV 

BAY  MINI-RANGER  NETWORK 


Figure  13.  SIAIION  COMBINATION  TREVOR-WATS-IIAVS-HANK 


58 


Figure  15.  STATION  COMBINATION  WATS-PACKARD-IIAVS-IIANK 
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